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I.  SUHMARY 

/ 

— ^  The  objective  of  this  work  is  to  develop  a  general  theoretical  framework  for 
calculating  Quctuations  of  signals  on  waves  propagated  through  random  media 
(WPRM)  and  to  apply  this  framework  to  sound  through  the  ocean;  light  through 
the  atmosphere;  radio  waves  through  the  ionosphere,  solar  wind,  or  interstellar 
plasma;  and  any  other  similar  case  of  waves  propagating  through  continuous 
media.  Comparison  with  real  data  is  an  important  aspect  of  the  effort. 

The  two  most  common  signals  sent  on  a  carrier  are  the  phase  and  ampli¬ 
tude  of  a  nearly  monochromatic  wave.  If  enough  bandwidth  is  available,  one  can 
send  a  pulse,  and  one  can  speak  of  the  intensity  and  arrival  time  of  that  pulse. 
The  technical  problem  is  then  to  explain  the  statistical  behavior  of  the  intensity 
and  arrival  time  in  terms  of  medium  fluctuations,  where  the  medium  is 
described  statistically,  usually  by  a  power  spectrum  covering  a  large  d)maniic 


range  of  scales. 


P  ^ 


The  eventual  practical  applications  of  an  understanding  of  YfPRM  to  science 
and  to  the  defense  department  are  myriad.  The  phase  of  a  light  wave  from  an 
astronomical  object  or  a  satellite  is  used  by  a  telescope  to  focus  to  a  detector; 
the  quality  of  the  focus  depends  on  the  state  of  the  atmosphere  -  a  random 
medium.  A  ground-based  laser  with  large  optics  attempts  to  focus  on  a  small 
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spot  for  a  period  of  time;  the  atmosphere  spreads  out  the  srot  by  its  action  on 
the  phase  of  the  wave.  Determination  of  pulsar  parameters  depends  on  observa¬ 
tions  of  radio  pulses  through  a  distorting  random  medium-interstellar  plasma. 
Communication  with  spacecraft  by  radio  pulses  using  radio  telescopes  depends 
on  coding  schemes  and  antenna  control  that  must  contend  with  effects  due  to 
the  solar  wind  and  the  earth's  ionosphere.  Communication  with  earth  satellites 
must  contend  with  the  ionosphere,  and  the  effects  of  a  disturbed  ionosphere 
(due,  e.g.,  to  nuclear  explosions)  must  be  predicted.  Probing  oceem  processes, 
from  large-scales  (Gulf  stream)  on  down  (internal  waves  and  microstructure) 
depends  on  understanding  acoustic  propagation  through  a  random  medium;  for 
example,  a  favored  method  is  to  send  pulses  over  long  range  and  observe  their 
arrivad-time  variations.  Detection  of  submarines  by  passive  acoustics  is  limited 
by  ocean  fluctuations  that  control  the  maximum  antenna  size  and  integration 
time  that  can  be  used  for  coherent  signal  Integration.  Active  sonar  for  commun¬ 
ication  has  a  limited  bandwidth  due  to  ocean  fluctuations.  Determination  of  the 
characteristics  of  earthquakes,  thought  to  be  important  for  earthquake  predic¬ 
tion,  is  done  solely  by  seismic  detection.  The  earth  through  which  the  seismic 
waves  travel  has  randonmess,  and  this  limits  the  information  that  can  be 
gleaned  from  seismic  signals.  The  same  effect  limits  our  ability  to  distinguish 
between  underground  nuclear  explosions  and  earthquakes,  and  thus  affects  our 
political  stance  vis-a-vis  nuclear  test  bans.  On  the  other  hand,  the  distortions  of 
seismic  waves  due  to  the  earth’s  random  properties  can  be  used  as  a  probe  of 
those  properties  and  hence  can  lead  to  a  better  understanding  of  earth  struc¬ 
ture. 

The  path-integral  method  for  treating  wave  propagation  has  been  success¬ 
fully  used  by  the  principal  investigator  for  the  analysis  of  many  experiments  in 
ocean  acoustics.'  This  method  is  therefore  utilized  extensively  in  our  work. 
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Another  school  of  wave  propagation  theory,  which  began  in  the  early  1960' s  in 
order  to  explain  experiments  in  light  transmission  through  the  turbulent  atmo¬ 
sphere,  utilizes  a  different  technique  involving  partial  differential  equations  for 
the  moments  of  the  wave  field.  The  moment  equations  and  the  path  integral 
have  now  been  used  by  enough  researchers  that  the  value  of  both  approaches  is 
appreciated,  but  the  relations  between  the  two  methods  has  remained  confusing 
to  many  people. 

A  first  step  in  YTPRM  is  characterizing  the  statistics  of  the  random  medium. 
In  this  report  we  will  speak  of  the  spectrum  of  medium  fluctuations,  which  would 
be  obtained  by  dragging  an  index-of-refraction  sensor  through  the  medium  and 
taking  a  Fourier  transform  of  the  resulting  time  (=  space)  series.  This  spectrum 
is  characterized  by  a  power  law  (e.g.  -5/3  for  Kolmogorov  turbulence)  that 
implies  much  more  variation  at  large  scales  than  at  small  scales.  It  is  likely  that 
a  sensor  dragged  in  different  directions  will  observe  different  spectra  even  on 
the  average.  In  that  case  we  speak  of  an  anisotropic  spectrum.  For  example, 
the  strength  of  the  ocean  spectrum  is  much  higher  in  the  vertical  than  the  hor¬ 
izontal  for  the  seune  wave-number.  Finally  it  is  important  to  know  that  the  spec¬ 
trum  cannot  continue  indefinitely  at  either  large  or  small  scales.  At  the  "outer 
scale”  the  spectrum  cuts  off  due  usually  to  finite  container  size  —  the  height  of 
the  atmosphere  or  the  depth  of  the  ocean.  At  the  "inner  scale”  the  spectrum 
cuts  off  due  to  physical  processes;  for  example,  viscosity  becomes  important  at 
scales  of  order  a  few  millimeters  in  the  atmosphere  and  ocean. 

This  technical  report  covers  the  third  six  months  of  our  contract  effort. 
The  next  several  paragraphs  summarize  the  technical  results  we  have  obtained: 
more  detadls  are  given  in  following  secUons.  The  work  has  been  carried  out 
under  the  direction  of  Dr.  Stanley  Flattd,  and  involves  effort  by  Dr.  Flattd,  senior 
scientist  Dr.  Frank  Henyey,  two  post-doctoral  researchers  Drs.  Dennis  Creamer 
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and  Rod  Frehlich,  and  a  graduate  student  (in  the  UCSD  Electrical  Engineering 
and  Computer  Science  Department)  Johanan  Codona. 

Our  progress  in  understanding  the  travel  time  of  pulses  in  random  media 
has  resulted  in  a  paper  ACCEPTED  for  publication.^  as  well  as  results  that  will 
lead  to  later  publication.  Pulses  sent  through  a  fluctuating  medium  arrive  ear¬ 
lier  or  later  than  they  would  in  the  absence  of  fluctuations,  depending  on  the 
particular  realization  of  the  medium.  The  variance  of  arrival  time  can  be  calcu¬ 
lated  by  straightforward  methods  in  the  geometrical  optics  limit  Our  dramatic 
new  result  is  for  the  average  arrival  time,  which  we  And  advanced  in  weak  media. 
Heretofore,  researchers  were  of  the  opinion  that  pulses  were  delayed  on  the 
average.  This  effect  is  of  little  importance  in  communication  applications  where 
the  average  arrival  time  is  usually  less  important  than  the  variance.  However, 
the  effect  can  be  important  in  probing  a  random  medium  for  large-scale  varia¬ 
tions  by  their  effect  on  average  travel  time.  Our  result  implies  a  possible  confu¬ 
sion  between  a  changing  turbulence  level  and  a  change  in  the  average  index  of 
refraction  on  a  large  scale.  For  example,  ocean  acoustic  tomography  attempts 
to  measure  the  warming  of  a  100-km-square  area  of  the  ocean  by  an  expected 
change  in  travel  time  of  about  20  ms.  However  we  And  a  cbeuige  in  average 
travel  time  of  about  10  ms.  due  to  an  internal- wavefield  that  has  no  average 
warming  at  all.  Ve  have  recently  studied  the  range  dependence  of  this  effect, 
and  have  found  that  it  grows  as  the  square  of  the  range.  This  implies  that  exper¬ 
iments  being  planned  in  the  1000-4000  km  region  will  have  major  difficulties 
sorting  out  the  effects  of  internal  waves  from  the  effects  of  large-scale  struc¬ 
ture.  Most  importantly,  the  determinations  of  internal-wave  effects  will  NOT  be 
contaminated  by  the  large-scale  effects. 

Moving  from  arrival  time,  which  is  related  to  phase,  we  discuss  amplitude  or 
intensity.  We  have  developed  a  method  for  calculating  the  spatial  correlation 


function  of  intensity  on  a  transverse  plane  through  a  receiver.  This  is  a  long¬ 
standing  problem  that  is  of  great  importance  in  using  wave  propagation  for 
probing  the  structure  of  a  random  medium,  because  measuring  intensity  is 
often  the  observation  that  can  be  made  most  easily.  In  addition,  an  amplitude- 
modulated  signal  will  be  degraded  by  intensity  fluctuations  due  to  the  medium. 
The  standard  theory  develops  a  series  solution  for  the  intensity  spatial  spec¬ 
trum.  The  first  few  terms  are  an  accurate  representation  of  the  small-wave¬ 
number  end  of  the  spectrum.  In  order  to  calculate  the  high-wave-number 
region  many  terms  of  the  series  had  to  be  evsduated.  We  have  determined  a 
dififerent  series  expansion,  whose  first  few  terms  give  the  high-wave-number  sec¬ 
tion  of  the  spectrum.  Hence  the  evaluation  of  the  full  spectrum  is  simplified 

considerably.  A  paper  describing  our  results  has  been  SUBMlTlKb  to  Radio  Sci- 
s 

ence. 

We  have  also  made  considerable  progress  toward  evaluating  the  intensity 
spectrum  for  an  arbitrary  source  distribution,  going  beyond  the  standard  pro¬ 
cedure  of  considering  the  special  cases  of  a  point  source  or  an  incident  plane 
wave.  Our  general  case  will  include  a  source  that  is  extended  over  a  large  aper- 

A 

ture.  We  also  treat  the  two-frequency  case.  An  example  of  a  coherent  source  of 
large  aperture  would  be  a  large-aperture  laser  beam  An  example  of  an 
incoherent  source  is  a  planet,  or  an  illuminated  satellite,  or  an  infrared  plume 
from  an  ascending  booster. 

Because  researchers  favoring  the  moment-equation  method  or  the  path- 
integral  method  typically  knew  only  one  of  the  methods  in  any  depth,  the  rela¬ 
tion  between  the  two  methods  has  been  a  mystery  to  many.  We  have  expended 
quite  a  bit  of  effort  to  understand  this  relation.  We  have  a  paper  in  REVIEW  by 
the  Journal  of  Mathematical  Ph}rsics  that  will  show  that  the  two  methods  are 
mathematically  equivalent,  in  much  the  same  way  that  the  Heisenberg  and 


Schr5dinger  approaches  to  quantum  mechanics  were  shown  to  be  mathemati¬ 
cally  equivalent.  A  better  analogy  for  those  familiar  with  quantum  mechanics  is 
the  equivalence  of  the  Schrodinger  and  Feynman  approaches  to  quantum 
mechanics.  The  equivalence  extends  to  the  equivalence  term  by  term  of  the 
series  solutions  for  the  intensity  spectra  mentioned  earlier. 

In  nearly  all  cases,  in  order  to  compare  theory  to  experimental  data  in 
WPRM,  we  must  use  a  model  spectrum  for  the  medium  fluctuations.  We  have 
developed  phenomenological  spectra,  as  a  function  of  wave  vector,  that  allow  for 

Q 

an  anisotropic  componf^nt  added  to  a  turbulent  isotropic  component.  This 
model  is  meaningful  both  for  the  ocean,  where  the  anisotropic  component 
represents  internal  waves,  and  the  ionosphere,  where  the  anisotropy  is  due  to 
electrons  preferentially  moving  along  magnetic  fleld  lines.  We  are  in  the  process 
of  calculating  intensity  spectra  in  the  weak  fluctuation  region  using  these  model 
spectra.  We  have  data  from  an  ocean-acoustic  experiment  that  will  be  used  for 
comparison  purposes;  the  experiment  utilized  10-70  kHz  sound  over  several  hun- 

A 

dred  meters  under  the  Arctic  ice.  We  should  note  that  the  weak-fluctuation 
regime  is  one  in  which  the  intensity  series  solution  for  the  low-wave-number 
regime  is  the  only  relevant  one. 

For  more  than  one  hundred  years,  eclipse  observers  have  noted  ''mysteri¬ 
ous”  bands  of  shadows  moving  on  the  ground  just  before  and  after  an  eclipse. 
Many  exotic  theories  of  these  shadow  bands  have  been  put  forward,  but  most 
observers  agree  that  they  are  probably  due  to  atmospheric  scintillation  that 
becomes  visible  when  the  crescent  of  the  eclipsed  moon  becomes  thin  enough  (a 
few  minutes  before  and  after  an  eclipse).  Johanan  Codona,  the  graduate  student 
associated  with  this  project,  has  made  the  first  systematic  application  of  WPRM 
theory  to  eclipse  shadow-band  observations.*®  He  explains  the  orientation  and 


contrast  of  the  bands  as  a  function  of  time,  and  describes  the  effects  of  eclipse 
geometry  and  the  importance  of  wind  direction.  He  relates  shadow-band  obser¬ 
vations  to  stellar-scintillation  observations.  One  important  conclusion  he  draws 
is  that  as  the  illuminated  crescent  gets  thinner,  the  shadow-band  observations 
probe  higher  into  the  atmosphere.  Recently  published  data^^  from  the  eclipse 
of  February  16,  1980  in  India  agrees  with  Codona’s  predictions.  Further  data 
from  the  annular  eclipse  of  May  30,  1984  in  Georgia  should  soon  be  forthcoming. 

We  have  begun  the  analysis  of  seismic  data  from  the  Center  for  Seismic  Stu¬ 
dies.  Two  nuclear  explosions  in  the  Soviet  Union  with  good  detections  on  the 
NORESS  array,  which  has  about  twenty  elements  spaced  out  to  a  few  kilometers 
have  been  obtained.  The  first  look  shows  rather  small  travel-time  fiuctuations. 
somewhat  at  odds  with  the  large  amplitude  fluctuations  that  have  been  sug¬ 
gested  previously.  We  have  looked  at  the  data  as  a  function  of  frequency  up  to 
about  20  Hz.  and  have  seen  no  obvious  systematic  differences  in  travel  time 
between  the  different  frequencies,  except  for  an  unusual  change  in  the  arrival 
structure  between  5  and  10  Hz. 

As  discussed  in  our  last  six-month  report,  this  period  and  the  next  few 
months  are  primarily  periods  of  preparation  of  manuscripts  for  publication.  We 
wiU  be  presenting  this  work  at  a  number  of  conferences,  including  a  meeting  of 
the  Union  of  Radio  Scientists  at  Vancouver  (June  1985). 

We  will  continue  to  be  directly  involved  in  comparison  between  theory  and 
experiment  over  the  next  year.  Spacecraft  data  from  the  Stanford  Center  for 
Radar  Astronomy,  satellite  data  from  SRI-International,  and  seismic  data  from 
the  Center  for  Seismic  Studies  (Virginia)  will  be  included,  along  with  our  ongoing 


ocean-data  effort. 


Our  longer-term  goals  will  include  the  implementation  of  our  new  theoreti¬ 
cal  results  into  computer  codes  for  calculation  of  general  phase  and  amplitude 
fluctuations.  Two  directions  arc  contemplated  that  will  require  large-scale  com¬ 
puting.  The  first  is  propagating  waves  through  individual  realizations  of  random 
media  to  compare  with  our  theoretical  results  and  to  extrapolate  those  results 
into  parameter  regimes  in  which  the  theory  is  not  valid.  This  propagation  can 
be  done  via  a  parabolic  wave  equation,  so  that  it  is  a  marching  solution.  The 
second  involves  evaluating  the  theoretical  formulas  which  involve  either  multidi¬ 
mensional  ordinary  integrals,  or  in  some  cases,  path  integrals.  We  have  begun 
some  computer  work  on  simulation  of  WPRM  using  our  VAX,  and  we  plan  to 
implement  the  code  on  an  IBM-PC  that  has  two  FFT  hardware  boards  that  should 
allow  uninterrupted  calculations  at  about  twice  the  speed  of  a  VAX.  These  efforts 
are  in  preparation  for  proposals  to  do  simulation  calculations  on  a  CRAY. 
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n.  EQUIVAlfNCE  OF  PATH  INTEGRAL  AND  MOMENT  BQUATIQNS 


Two  approaches  have  yielded  many  results  for  WPRM.  These  are  the 
moment  equation  approach  and  the  path  integral  approach.  Although  they  both 
start  with  the  parabolic  wave  equation,  approximation  are  made  in  each  to 
derive  the  fundamental  equation  or  expression  for  any  statistical  moment  of  the 
wave  field.  These  approximations  are  in  the  nature  of  assuming  that  correlation 
lengths  are  not  too  long.  The  relationship  between  the  two  approaches  has  not 
been  clear. 

We  show  that  the  two  approaches  are,  m  fact,  equivalent.  By  using  the 
method  that  Feynman  used  to  show  his  path  integral  was  eqmvalent  to 
Schrodinger's  equation,  we  derive  an  equation  equivalent  to  the  path  integral, 
and  find  that  this  equation  is  the  moment  equation  used  in  practice  as  a  starting 
point  in  WPRM  studies. 

In  some  applications,  both  approaches  yield  expansions  of  the  moment  in 
some  vaiiable.  Since  the  approaches  are  equivalent,  the  expansions  are  equal 
term  by  term.  In  other  applications,  one  or  the  other  approach  might  yield  a 
desired  expansion  or  approximation  more  readily  than  the  other.  For  example, 
corrections  to  the  "semiclassical”  result  are  more  easily  obtained  from  the  path 
integral. 

The  following  manuscript  is  our  finished  paper  on  this  subject  It  has  been 
submitted  for  publication  in  the  Journal  of  Matbematiced  Physics. 
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L  IntrodiaeUco 


Many  problems  in  wave  propagation  through  random  media  concern 
phenomena  in  which  there  is  no  signiAcant  backscalter,  m  that  a  parabolic 
approximation  may  be  made  to  the  wave  equation.^'1  In  these  cases  a  further 
approximation,  called  the  Markov  approximation. leads  to  relatively  tractable 
mathematical  expresnons  for  momenta  of  the  held  that  can  be  used  for  practi¬ 
cal  calculatims.  ^ro  quite  different  f<x‘malisms  have  been  used  in  this  context: 
the  moment-equation  and  path-integral  techniques. 

A  path-integral  expression  for  a  general  moment  of  the  Aeld  of  a  wave  pro¬ 
pagating  through  an  inhomogeneous,  amsotropic  medium  in  the  presence  of  a 
deterministic  background  refractive  index  has  been  derived.I*l  and  the  eiq>res- 
sicxi  has  been  used  for  speciAc  calculations.^^*-*^ 

Moment  equations  in  coordinate  representatimi  have  been  derived  for 
homogeneous  isotropic  media  in  the  absence  ot  a  detemunistic  background. 
Treatments  of  inhomogeneity,  anisotropy,  and  deterministic  background  by 
moment-equation  techniques  have  herettrfore  been  conAned  to  special  cases 
involving  the  ffrst  and  second  momenta. ’ 

We  present  here  general  moment  equations  in  coordinate  representation 
that  account  for  inhomogeneity,  anisotropy,  and  deterministic  background,  but 
require  the  Markov  approximatioa  We  derive  these  equations  using  the  time- 
ordered-product  method  of  Van  Kampen.^*!  which  also  provides  a  derivation  of 
equations  that  are  valid  under  conditions  more  general  than  the  Markov  approxi¬ 
mation.  The  modiffed  equations  are  more  complicated  than  those  that  require 
the  Markov  approximation:  a  special  case  was  previously  derived  by  Besieris  and 
TapperLf*®^ 

We  also  show  that  our  new  general  moment  equations  derived  under  the 
Markov  approximation  are  mathematically  equivalent  to  the  path-integral 
expressions  for  the  moments  that  have  been  previously  presented.  Thus,  the 
two  popular  formalisms,  under  the  Markov  approximation,  are  not  different  in 
content 


The  plan  of  the  paper  is  as  follows:  in  Section  D  we  establish  notation, 
present  our  new  moment  equations,  and  present  path^tegral  expressions  for 
the  moments  in  similar  notation.  In  Section  SI  we  estabUsb  the  mathematical 
equivalence  between  the  two  techniques.  In  Section  IV  we  present  the  derivation 
of  our  moment  equations,  and,  along  the  way.  derive  the  modified  equations.  In 
Section  V,  for  completeness,  we  rederive  the  path*integral  eq>ressions  for  the 
moments.  In  Section  VI  we  comment  on  the  use  of  different  coordinate  systems 


0.  NotaUco  and  MarkofnApprMdmftUoD  RMiilta 

Consider  waves  IreveUing  predominantly  in  the  s  direction.  Let  f  be  e 
transverse  coordinate  (e.g.  two-dimensional,  but  in  fact  general),  and  b  be  a 
reference  wave  number  {k  =  2no/  Cq).  where  u  is  the  wave  frequency  and  Cq  is 
a  reference  wave  speed).  Express  the  full  wave  field  as 

u(f.s.t)  =  f(d.s)exp[ib(s  -Cof))  (1) 

Let  the  wave  speed  (a  function  of  position  only)  be 

C(f.s)  =  Ce[l-2£/o(f)-?M(*.e))-»  *  Co[l  ♦  t;,(f)  ♦  (2) 

where  C/q  represents  the  detenniniaUc  background  and  n  represents  the 
fluctuating  random  medium,  assumed  to  be  a  realisation  of  a  zero-mean  Gaus¬ 
sian  process. 

Then,  the  parabolic  equation  (in  rectufigular  coordinates)  for  the  reduced 
wave  function  ^  is : 

tte.f  =  -  ♦  k^fi(*.t)f  (3) 

where  V*  is  the  transverse  Laplacian. 

A  moment  T  is  the  ensemble  ezpeciaticm  value  of  a  product  of  i^s  and  ^*’s 
where  each  f  or  is  evaluated  at  a  different  position  tf  and  wavenumber  k^. 
We  write,  in  abbreviated  form. 

~  ^1^1  '  ■  ■  i^m*n  ^  (4) 

Define  an  operator  Lq  such  that 

The  terms  that  apply  to  the  f  s  use  the  phis  sign  and  those  that  apply  to  the 
use  the  minus  siga  The  subscript  i  requires  that  Vf  operate  only  on  if  and 

Define  the  important  cnnbination  of  fluctuation  quantities  as 
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H{9)  *  y!" (6) 

Our  general  moment  equation  under  the  Markov  approximation  can  be 


written 

-  X  /*■<*(« C) 

where  )  is  obtained  by  evaluating  itf  (e  )  with  all  the  at  s  shifted  by  the 

transverse  distance  that  a  determimstic  ray  through  moves  in  travelling 

from  c  to  s*  (see  Figure  1).  In  other  words  if^|(s‘)  is  evaluated  at  point  B:  Le. 

=  W)  where  the  ray  is  forced  to  go  through  f^(s).  Ibe  particular  ray  is 
determined  not  Mily  by  the  local  position  (^j.s).  but  also  by  the  initial  ccxidi- 
tions  on  the  moment;  for  example,  the  location  of  a  point  source,  or  the  direo- 
tion  of  a  plane  wave.  The  unphysical  assumption  of  d<!lta-corr elated  medium 
fluctuations  along  the  propagation  direction  would  imply  that  (s' )  would  be 
evaluated  at  point  C:  Le.  f/(s)  (and  s').  In  the  isotropic  case  (or  in  the  case  of 
propagation  altxig  a  principal  axis  of  the  anisotropy)  the  difference  between 
evaluating  ^/(*)  >*  negligible,  and  the  delta<correlated 

assumption  is  adequate.  In  the  anisotropic  case,  the  necessity  of  defining  the 
unperturbed  ray  makes  (7)  somewhat  complicated  to  apply  for  general  initial 
conditions.  However,  since  (7)  is  a  linear  equation  superposition  can  be  used 
whether  the  source  is  a  point,  an  incident  plane  wave,  or  an  arbitrary  coherent 
or  incoherent  sum  of  point  sources.  Equation  (7).  which  is  one  of  the  principal 
results  of  this  paper,  is  derived  in  Section  IV. 

We  now  turn  to  the  path  integral  method.  E(pjation  (3)  has  the  formal  solu> 


y  =  /Pf(s)e«  (8) 

where  f  means  integration  over  paths,  i{t)  is  a  transverse  vector  indi¬ 

cating  the  position  of  the  path  at  s.  and 


5  =  *  j  -  t7o(*)  -/*(f.  «) 


In  order  to  obtain  a  fiven  moment,  expresiiona  like  (B)  (or  its  complex  eonju* 
gate)  are  multiplied  together,  and  the  ensemble  average  is  taken: 


^  f  >  (10) 

Tbe  Markov  approximation  yields  (See  Section  V): 

(11) 

We  show  in  the  next  section  that  the  moment  equations  (7)  and  the  path  integral 
expressions  (11)  are  mathematically  equivalent. 
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ID.  Equlralciice  of  Paih  hiUfral  and  Moamit  Cquatlou  ondar  Um  MaHkor 
ApproilmaUeo 


We  follow  the  technique  that  Feynman^''^  used  to  show  that  his  path«integral 
expression  for  nonrelalivistic  quantum  mechanics  is  equivalent  to  the 
Schrfidinger  equation.  Ilie  key  to  this  demonstration  is  an  understanding  of  how 
the  important  paths  behave  transversely  as  they  move  in  *  from  a  particular 
point.  Feynman  found  that  these  paths  resembled  random  walks  in  that 

|f(s')-f<s)|-(s'-s)^  (12) 

as  s'  gets  close  to  s.  Given  this  beha^or.  it  is  easy  to  expand  (11)  in  a  Tajdor 
series  and  obtain  a  dilTerential  equation  which  will  turn  out  to  be  (7).  We  give 
the  demonstration  of  (12)  in  the  Appendix. 

The  path  integral  is  defined  as  the  limit  of  an  integration  over  a  set  of 
"phase  screens."  These  screens  are  at  values  s  Nit.  The  derivative  ^"nt 

z  s  tft  is  defined  as  it{zti  +  ds)  -  f{xfi))/  it  ~  The  limit  ds-*0  is  taken 

after  the  integrals  are  evaluated.  The  differential  equation  is  obtained  by  com 
sidering  the  integral  over  the  very  last  phase  screen.  The  last  integral  in  (11) 
can  be  written  in  terms  ot  I  j'  =  I  f{F  -  it)  and  f  f  -  i  j{F).  Also,  we  define 
6ij  -  tj Then  F,,,  can  be  expressed  as: 


(13) 


v^{\r\,R-iz) 


where  denotes  the  set  of  m-fn  f^'s.  The  first  term  in  the  exponent, 

iJtf  I6ii)* 

— .  is  0(1)  for  small  it,  because  of  (12).  The  exponent  of  the  remain- 

6  OS 

ing  terms  can  be  expanded,  since  they  have  an  explicit  ds,  as  well  as  higher 
order  terms.  This  results  in 


4-»-« 


(14) 


I 


ft 


exp  I  tifcy  (<!/)•/ 2«ti) 


1-d.  (£±0:^  C/«(y?)+  |-7  d.'  <ll(/?)i#^,(t  )» 

i  “• 


U-<«)  ♦  o(d.^) 


We  now  have  a  relationship  between  the  moment  at  R  and  the  moment  at 
R-i*.  which  we  derived  from  our  path-integral  expression.  But  since  the 
moment  is  a  differentiable  function  we  can  find  another  relationship  by  Taylor 
expansion  as  follows: 


1  -  ds  dj,  -  £df,.V^  + 


(15) 


♦  0(ds*) 


Substituting  (15)  into  (14)  we  find 


■^s  d/f—iz 


^±ikfU^iR)*^fjb’  <  if(/?)ir^,(s’)  > 


2. 

r«,(j*U)  +  o(ds2) 

The  term  linear  in  is  odd  in  6if  and  therefore  gives  zero  due  to  the  6lf 

integral.  The  term  that  is  quadratic  in  iff  can  be  integrated  by  parts,  yielding, 
to  order  ds: 


•  1  ♦  4.  - »,  -  k;  -^vf  -  z  *«,  wo,(*) 

C  ^  tWj  ^ 

-|-/(i.‘<i/{^)i/^,(«)>|r«*«fU)  (17) 

The  only  way  (17)  can  be  true  for  all  ds  is  for  the  coefficient  of  6t  within 
the  curly  brackets  operating  on  Pmi  to  give  zero.  Therefore,  setting  R-**  • 

(»B) 

which  is  identical  to  (7).  as  required.  Thus,  we  have  derived  the  moment  equa¬ 
tion  (7)  from  the  path  'mtegral  expression  (11).  This  shows  that  the  path- 
integral  expression  (11)  is  a  solution  of  the  moment  equation  (7)  and  hence  the 
two  techniques  are  equivalent 
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IV.  lloai«ii*EquaUoo  Dorifatloo 

We  derive  our  moment  equetions  by  the  method  of  Van  Kampen.^*)  The 
advantage  of  his  method  is  that  the  physical  basis  for  each  approximation  is 
readily  apparent  He  bases  his  method  on  techniques  that  were  developed  for 
quantum  mechanics. 

We  shall  And  that  the  Markov  approximation  requires  that  the  dimensionless 
number  Qlif  be  small  where  is  the  medium  correlation  length  in  the  direc¬ 
tion  of  the  wave  propagation,  and  if|  is  the  “typical"  value  of  if.  deAned  by  (6) 
and  called  the  “interaction  strength."  For  the  Arst  moment  Jl  -  kfi,  but  for 
higher  moments  if  is  the  sum  and  diAerence  of  a  number  of  fc/i's  at  different 
positions,  and  with  different  values  of  k. 

Ve  start  with  the  parabolic  wave  equation  (3)  and  the  deffniticm  ot  L9  and  if. 
and  write: 

Va  •  =  (io  ♦  if)Vi* Vs  •  •  V*-  (19) 

The  “interaction  representation"  is  deAned  by: 

(ViW  •  Vm>ii)/  -  VrVs*  •  (20) 

and 

if/(s)  =  (21) 

With  these  deAnitions,  (19)  becomes 

i9.(VrVs*  V»»i.)/  =  fm*n)l  (22) 

This  equation  is  linear  and  has  the  formal  solution 

(Vr  Vs'  Vmtu)/  =  Texp  {-if  ’)dz  *)  T^{0)  (23) 

0 

rinn(O)  is  the  initial  conditioa  The  "time-ordering"  symbol  T  requires  explana¬ 
tion.  One  notices  that  if/  is  an  operator,  not  just  a  function  of  space.  if/(S|) 
and  Mi(tz)  do  not,  in  general,  commute.  If  they  did  the  solution  of  (22)  would  be 
given  by  (23)  without  the  T  symbol  The  T  symbol  means  that  a  product  of 
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operators  to  the  rifbt  is  not  ap^ied  in  the  usual  order,  but  in  such  a  way  that 
operators  with  smaller  values  of  e*  are  to  be  applied  flrst  Thus  there  i»  an  ord¬ 
ering  in  s.  (The  T-symbol  was  invented  for  solving  problems  in  quantum 
mechanics  where  the  analog  of  the  longitudinal  direction  e  is  the  time.)  For 
example. 

rerp(-i/lf/(e)(to)  =  (24) 

0 

[  Texp( -i  /  if/(*’ )d*' )  ]  [  7*exp(-i  /  j[//{e’ )d*’ )  ] 

•i  « 

Another  example  is 

T  Jf  J//<s')de']*  =  (-i)*  Uiitt)  lf;(e,)d*,  de*  .  .  .  da,  (25) 

where  the  integration  region  on  the  right  side  of  (25)  isO<Cj  <  .  .  .  <  <  x, 

which  is  k\  times  smaller  than  that  of  the  left  side,  cancelling  the  factor  of  k\. 
Using  either  (24)  or  (25).  one  readily  checks  that  (23)  is  a  formal  solution  of 
(22). 

We  are  assviming  that  il  is  a  Gaussian  process.  The  result  that  the  expecta¬ 
tion  of  the  exponential  of  a  zero-mean  Gaussian  random  variable  is  the  exponen¬ 
tial  of  half  the  variance  follows  from  combinatorial  factors  and  remains  true  tar 
a  time-ordered  exponential.  Thus 

(r^)/  =  Texp  -  |-<  if/(*’)d*'  >  r^{0)  (26) 

Although  this  is  a  formal  expression  for  Paw  it  ix  not  immediately  useful  for  cal¬ 
culations,  since  there  is  no  simple  algorithm  for  evaluating  a  time-ordered 
exponential  (in  contrast  to  a  normal  exponential).  Van  Kanq>en  proceeds  by 
differentiating  (28): 

a.(r-«)/  =  (27) 

-  r  <  Mi(x)  Jf  dM'Uiix’)  >  exp  -  [j  i/;(s")ds"  >  r^{0) 


a-a-ia 


The  itf/(e )  het  the  lergeii  c .  so  it  is  written  in  the  proper  ordered  position.  The 
if/(e' )  that  it  is  eorreleted  with,  however,  might  occur  anywhere  relative  to  the 
if/(s**  )*s  in  the  exponenUaL  If  L^ilf  oc  i.  very  little  error  is  made  by  assuming 
Uiat  the  first  two  if/'s  are  in  the  proper  order,  so  that  the  T  symbol  can  be 
brought  through  the  first  expectation  value,  yielding: 

a.(r««)/  =  -  <i//(e)|dsjf;(s)>(r^),  (28) 

This  may  be  shown  by  expanding  the  exponential  operators  in  (26)  or  (27) 
and  discussing  the  order  of  M‘s  in  each  term.  The  Nth  term  in  the  expansion  has 
2N  occurrences  ct  lii,  and  is  of  a  magnitude 

<{f  UidM  Y>"/  m  (29) 

0 

where  typical  eigenvalues  of  the  operators  are  implied.  The  terms  beyond 

N»»4<ifMidz)F>  (30) 

0 

become  negligible  compared  to  the  original  exponential,  so  we  have  to  deal  with 
at  most  N  pws  of  i#/'s  from  source  to  range  s .  The  two  lf/*s  in  a  correlated  pair 
must  be  within  Lp  of  each  other  to  give  a  nonmro  cmrelation.  The  number  of 
pairs  may  be  estimated  as 

(31) 

where  is  a  typical  value  of  if/.  (See  Figure  2  for  a  schematic  representa¬ 
tion.)  Our  approximation  reduces  to  sajring  it  is  unlikely  to  find  a  third 
occurrence  of  an  if/  in  between  a  pair  that  are  within  of  each  other.  This  pro¬ 
bability  is  roughly 

Probability  w  (32) 

Thus  if  the  fiuctuations  are  weak  enough  (if|  small),  the  approximation  is  valid, 
and  (28)  is  justified. 
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We  call  (26)  “first  order  perturbation  theory.'*  In  typical  situattons.  t  is 
much  larger  than  Lf,  and  the  lower  limit  can  be  replaced  by  »  «•.  making  the 
equation  independent  of  the  source  position.  Moreover,  the  integral  from  • 
to  s  can  be  replaced  by  half  the  integral  from  -  to  when  the  correlation  is 

a  much  slower  function  of  ^s  -fc')  than  of  s  The  result  is  used,  not  in  the 

interaction  representation,  but  in  the  original  representation.  The  exponentials 
of  (20)  and  (21)  are  removed,  giving 

d.r^i*)  *  (33) 

-  /  ds*  <  Jf(s)e“*^»‘'  lf(s*)e*^“  >  r,*(s) 


For  the  second  moment,  this  equation  is  related  to  an  expression  of  Besieris  and 
Tappert^'*’^  Although  their  work  was  for  the  second  moment,  we  can  generalize  it 
directly,  therefore  in  the  rest  of  our  comments  we  treat  the  general  moment 
Tiim  where  Besieiis  and  Tappert  treated  only  Pu.  Their  equation  3.2  was 
expressed  in  a  Fourier-transformed  domain,  but  can  be  expressed  in  our  nota¬ 
tion  as 

This  equation  is  equivalent  to  (33)  to  order  Lf  M*.  Unlike  (33).  (34)  implies  a 
“memory'*  effect  in  which  the  gradient  of  the  moment  depends  explicitly  cm  the 
moment  at  all  previous  s*s.  The  Markov  approximation  leads  to  (7).  which  elimr 
inates  the  memory  effect  and  requires  only  a  correlation  function  of  the  medium 
along  a  specified  (shilled)  direction.  Besieris  and  Tappert  pointed  out  that  a 
weaker  approximaticm,  called  the  “long-time  Markov"  approximation  leads  to  a 
local  (non-memory)  equation  (their  ecpiation  3.3),  that  in  our  notaticm  is 
expressed  as 

«.r«.  =  -1/ (35) 

We  are  only  ccmsidering  situations  in  which  the  parabolic  wave  equation  is  valid. 
It  has  been  shown  that  in  that  case  the  long-time  Markov  is  valid.^'*)  and 
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therefort  (35)  if  4i  vilid  as  (33). 

Because  Le  i*  an  operator,  the  integrals  in  (33>35)  involve  the  medium 
correlation  function  in  all  directions,  or.  in  the  Pourier*transform  domain, 
require  a  scattering  kernel  as  a  (unction  of  scattering  angle.  The  Markov 
approximation  to  (33)  consists  of  simpUtying  the  deterministic  propagation 
operator  e^^^*  for  s  -  s'  on  the  order  of  Lf.  Instead  of  correlating  if(s) 
with  all  possible  transverse  positions  of  Af(s'),  the  Markov  ai^roximation 
corresponds  to  choosing  only  one  transverse  position  for  if  (s').  (See  Figure  1. 
where  point  A  represents  an  arbitrary  transverse  position.)  If  the  wave 
represented  by  raM(>)  were  the  unperturbed  solution,  then  deterministic  pro* 
pagation  would  move  the  phase  in  the  direction  of  the  unperturbed  ray.  If  the 
wave  energy  is  traveling  close  to  the  unperturbed  ray  this  operator  retains  its 
behavior  to  first  approximatioa  As  a  result,  deterministic  propagation  approxi* 
mates  a  shift  along  the  unperturbed  ray  to  point  B.  Le..  f(s')  s  f,^(c‘).  where 
the  ray  is  forced  to  go  through  f(s).  Hence  e“*^***“^^if(s’)e*^*^*'^^  can  be 
approximated  by  Af«^/|(s‘).  Tlus  is  the  appropriate  definition  of  the  Markov 
approximation  (rather  than  assuming  the  medium  is  delta'correlated  along  the 
s  axis)  and  it  immediately  yields  (7)  from  (33).  In  practice,  instead  of.  using  the 
actual  unperturbed  ray.  the  tangent  to  the  ray  at  s  is  often  used. 

If  the  delta-correlated  assumptimi  were  made,  it  would  correspond  to 
evaluating  Af^i(x')  at  point  C,  which  is  strictly  valid  only  if  there  is  a  single 
unperturbed  ray  travelling  along  the  s-axis.  If  the  medium  fiuctuations  are  iso¬ 
tropic.  the  correlation  any  point  at  s'  with  the  point  P  at  s  will  give  the  same 
result  because  of  the  parabolic  approximation,  and  hence  the  delta-correlated 
assumption  is  as  good  as  any  other  choice.  However,  for  an  anisotropic  medium 
it  is  important  that  point  B  (and  hence  (7))  be  used,  even  when  the  Markov 
approximaUon  is  invtdced.  Note  that  (7)  can  be  used  in  the  presence  of  a  deter¬ 
ministic  background  refractive  index. 

The  difference  between  (33)  and  (7)  can  be  caused  by  directions  different 
from  the  unperturbed  ray  becoming  imjportant.  A  transverse  wavenumber  kf, 
coming,  for  example,  from  M,  causes  \ix  angle  to  change  by  dd  ^  kf  /  k.  A 
transverse  error  in  position  of  about  kfLp  /  k  'a  made  by  assuming  the 
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dirtciion  of  the  unperturbed  rey.  Thus,  in  order  (or  the  Markov  appronmation 
I  to  be  valid,  it  is  required  that  kjLf  /  k  dc  Ly.  where  is  the  transverse  scale  of 

\  concern.  Since  the  Markov  approximation  (ails  at  sufficiently  small 

k  tLf/  Lf,  The  parameter  a  ~  kLf/ Lf  introduced  by  Beran  and  McCoy^'*I  and 
discussed  further  in  Flattd^^l  reflects  these  considerations.  For  small  a.  one  can 
use  (33)  or  its  equivalent 


{ 


V.  PfeUHnUgnJ  Dtrtfalleo 

•  W«  rvcapiiulaU  the  dcrivetion  of  the  potb-integrel  exprceiion  (7)  from  (10). 
Using  the  essumed  Gaussian  behavior  of  the  fluctuations,  we  obtain  from  (10) 

r-  =  /  n  »/<•)•'  (“) 

•  /.I 

where  5q  is  the  part  of  5  in  (9)  that  does  not  involve  $t,  and 

V  -  -  ^f  da  dK‘  <  lf(e)lf(s')  >  (37) 

Ihe  eiqpression  (36)  is  an  exact  represeniaiion  of  the  moment  of  the  solution  of 
the  parabolic  equatitm  with  Gaussian  fluctuations.  It  is  not  used  in  practice  as  it 
^  stands  because  V  depends  on  the  paths  at  two  values  of  s,  namely  s  and  s' . 

The  Markov  approximation  for  the  path  integral  comes  from  assuming  that 
the  paths  do  not  stray  far  in  transverse  space  over  a  distance  Lp;  they  all  move 
approximately  parallel  to  the  unperturbed  ray.  Thus,  in  the  Markov  approxima* 

•  tion 
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K  r  -  dz  ds-  <  lf(s)jyuv,(s')  >  (38) 

which  only  requires  knowledge  of  the  path  at  t.  The  final  result  (11)  follows 
directly. 
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VI.  CoordliiaU  aystomi 

llomcni  equations  eon  bo  formuloiod  in  •  variety  of  coordinate  systems, 
while  path  integrals  require  a  rectangular  coordinate  system.  There  has  been  a 
fair  amount  of  effort  expended  on  using  polar  coordinate  systems,  especially  for 
point  source  problems. 

The  same  results  (for  point  sources  among  others)  can  be  obtained  in 
either  polar  or  rectangular  coordinates.  Thus,  the  results  of  Shishov^**^  on  the 
intensity  correlation,  derived  in  spherical  polar  coordinates,  can  be  seen  to  be 
identical  (after  an  appropriate  transformation)  to  the  results  of  Codona  ef 
derived  in  rectangular  coordinates.  It  was  necessary  for  Shishov  to  make  small 
angle  approximations  in  addition  to  the  parabcdic  approximation  of  dropping  the 
second  derivative  in  the  propagatim  direction,  whereas  Codona  ef  of.  only 
require  the  single  parabolic  apiM-oximaticm. 


W«  havt  derivtd  moment  equetione  in  coordinate  repreeenUUon  under  the 
Markov  approximation  that  apply  in  anisotropic.  inhomo(eneous  media  with 
deterministic  background.  The  derivation  shows  the  relationship  between  these 
moment  equations  and  modified  equations  that  are  vaiid  under  approsimaUons 
weaker  than  Markov,  the  second-moment  equation  of  Besieris  and  Tappert  is  a 
special  case  of  these  modified  equatimis. 

In  a  hierarchy  of  approximations  we  begin  with  the  parabolic  wave  equation 
itself.  A  path  integral  with  non-local  exponent  can  be  written  as  an  exact  solu¬ 
tion,  although  it  is  not  yet  useful  in  practice.  The  next  level  is  the  approKkna- 
tion  that  the  interaetimi  strength  over  a  cmrrelation  length  is  small-this  "first- 
order  perturbation  theory"  leads  to  the  modified  moment  equations,  and  in 
homogeneous,  isotropic  media,  to  the  standard  moment  equations  and  path- 
integral  expressions.  In  anisotropic,  inhomogeneous  media,  however,  a  hirther 
approximation  is  necessary  to  obtain  the  moment  equations  and  path  integral 
expressions.  This  further  approximatimi  is  that  the  significant  flow  of  wave 
energy,  or  the  important  paths,  are  parallel  to  the  unperturbed  ray;  we  call  this 
the  Markov  aiproximation  because  its  violatimi  implies  the  aiq>earance  of  corre¬ 
lations  between  successive  scatterings.  We  have  shown  that  the  moment  equa¬ 
tions  and  the  path-integral  expressions  for  the  moments  are  mathematically 
equivalent  under  the  Markov  approximation.  Thus  the  two  formalisms  have 
exactly  the  same  physical  content.  In  an  anisotropic  medium,  the  moment 
equation  involves  a  shift  (^ration  to  calculate  the  medium  correlation  function 
along  the  unperturbed  ray;  this  form  of  the  moment  equation  has  not  been  given 
before. 

We  have  also  pointed  out  that  aU  appropriate  formulae  can  be  derived  in  a 
rectangxilar  coordinate  system  (even  fw  point  sources). 


ildmowledgement  This  work  was  supported  by  DARPA. 
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Appendix 


We  must  show  that  the  scaling  |Afj  |  -  (ds )  ^  holds  for  integrals  of  the  form 


(^l) 


If  /*  is  expandable  in  a  power  series  (even  if  the  radius  of  convergence  is  zero) 
this  result  fcdlows  immediately.  One  expands  F  and  integrates  term  by  term. 

i. 

obtaining  a  power  aeries  in  (ds)^.  By  standard  methods  in  the  theory  of  asymp> 
totic  e]q>ansiona.  only  the  low  order  teims  need  to  be  retained  as  dt  -*0. 

For  singular  functicxis,  a  demonstration  is  not  as  simple  .  One  may  worry 
about  cancellations  between  terms  in  the  exponent,  since  the  signs  might  differ. 

We  will  content  ourselves  with  a  demonstration  in  the  ease  likely  to  arise  in 
practice.  It  is  common  to  model  a  random  medium  as  having  a  power  law  struc¬ 
ture  function.  Thus  as  two  x‘s  become  equal,  a  singularity  |Xi  -  \*  with  p  >0 

might  occur  in  the  integrand.  In  order  to  have  possible  cancellations  in  the 
exponent,  we  assume  that  kt-kf  s  k,  wd  the  exponential  factor  is 
exp(ik(dx(^  ~dx/)/2ds).  We  assume,  for  simplicity,  that  Xi  and  Xf  are  one¬ 
dimensional:  higher  dimensional  singularities  are  effectively  weaker. 

DeAne  v  -  (dx*  +  6zj)/Z.  ft  =  6xi  —  6Xj.  o  =  Xj  -x,-.  The  singularity 
hvm  the  previous  step,  x^'  =  Xf  -  iXf  is  |/a  ~  a  I'-  ‘Hie  integral  to  be  evaluated 
is 

f  dftdv  e**»^**  f(a./t,i/)  Im  “  «l'  (A.2) 


We  would  like  to  ignore  the  p  dependence  in  f.  However,  spurious  large-/z 
contributions  would  arise,  even  though  we  are  only  interested  in  contributions 
from  p  close  to  a.  To  drop  the  p  dependence  of  f  and  also  to  simplify  the 
analysis,  we  introduce  a  convergence  factor  exp(-a(/A*  4-  i^)/dx'  '*).  As  l«ng  as 
p,v  ~  dsA  this  factor  does  not  change  the  integral  as  dt  -»0  (we  are  assuming 
e>0).  Conversely,  if  the  integral  in  the  limit  dt -*0  does  not  depend  on  a  and  e. 
then  p  and  v  are  of  order  dxA 
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The  inUgrel  is  then 

/  *  fittdv  exp[-e(^  ♦  i^/dt  *"*]  \n  -  a|>  f(ci.i^)  (A.3) 

The  ti  integral  can  be  done: 

/  s  Ciy^di'e***'**  exp[-a(a*  ♦  (A,4) 

where  J/  is  a  confluent  hypergeometric  function  and  C}  is  a  constemt  indepen¬ 
dent  of  ds  .e.  and  a  (as  are  Cc  and  Cs.  below).  The  hypergeometric  functicm  has  a 
part  that  behaves  as  the  eoqionential  of  its  argument  for  large  (positive)  values 
of  its  argument,  a  part  that  falls  as  a  power  (since  (p<ff)/2  is  positive)  and  a 
part  at  small  values  of  the  argument.  These  last  two  parts  can  be  combined  into 
a  bounded  part.  We  show  that  the  exponential  part  gives  the  leading  behavior 
and  the  bounded  part  is  a  higher  power  of  dc. 

The  contribution  /|  from  the  exponential  asymptotic  part  ot  Jii  is 
I^  s  Ca/d./e*-'''**  exp[-a(o*  + u*)/d*‘-*l  (A.5) 

f(o.«/)  a-f»V2g,(i*p)(i-eV2  exp 
The  exponential  from  M  cancels  much  of  the  first  two  exponentials: 

which  can  be  done  explicitly.  Only  the  first  term  in  the  exponential  survives  as 
ds  -*0.  The  result  is  independent  of  a  and  e,  and  is 

Ix  -  Cad*f(a,0)|a|'  (A. 7) 

exactly  as  would  be  obtained  from  the  Taylm*  series  expansion  for  1. 

We  now  turn  to  the  contribution  /a  frmn  the  bounded  part  of  M.  We  show  /a 
has  a  higher  power  of  de  than  /|.  We  can  set  t  to  any  positive  value.  At  large  t 
we  depend  on  the  fact  that  averages  to  aero  for  for  any  positive 


(>»-f8idt*att) 

4ade>*« 
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0.  but  it  would  be  necessary  to  examiM  the  detailed  behavior  of  i/  to  use  this 
fact  On  the  other  band,  for  small  enough  s.  it  suffices  to  bound  the  integral  by 
the  integral  of  the  absolute  value  o^  the  integrand.  The  convergence  factor  pro¬ 
vides  a  cutoff  at  Thus  f  dv  t{a,v)  gives  a  contribu¬ 

tion  scaling  like  6s  Thus  It  is  bounded  by  an  expression  which  scales  as 


As  long  as  we  have  chosen  e  small  enough,  the  expiment  of  6s  is  larger  than  1, 
and  It  can  be  neglected  relative  to  1%.  Thus  we  have  established  the  necessary 
scaling  of  fi  and  v  even  in  the  singular  case. 
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Figure  lb. 


Figure  2a. 


Figure  2b. 


Moment*«qualion  axpreuion  of  the  Markov  approximaUoa  Ibe 
correlation  should  be  taken  between  a  point  at  t  (point  P)  and  an 
arbitrary  point  at  s’  (point  A).  Instead  it  is  taken  with  the  point  B. 
obtained  by  extrapolating  along  the  unperturbed  ray  from  F.  The 
assumption  of  delta*correlated  medium  fluctuations  leads  to  the 
incorrect  formulation  of  correlations  between  points  P  and  C.  The 
dashed  lines  indicate  the  idea  of  a  scattering  as  a  fiinction  of  angle 
from  point  P. 


Path-integral  expression  of  the  Markov  approximation.  Tbe  general 
path  at  s'  (point  A)  is  approximated  by  the  path  at  s  extrapolated 
altmg  the  unperturbed  ray  (p<tot  B). 


Typical  s  values  of  the  interactions  from  a  Taylor  series  term  in  (9) 
are  indicated  by  s’s.  Dashed  lines  show  which  interactions  are 
correlated.  It  is  assumed  that  Z^JV‘«1. 


A  portion  of  a  contribution  to  (9)  which  is  improperly  ordered  in 
"first  (x-der  perturbation  theory."  Such  centr&utions  are  small  il 
£/if,*«l. 
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ni.  INTDISmrCORRQATIONS 

The  following  section  is  a  paper  on  a  new  treatment  of  the  fourth  moment 
that  allows  calculations  in  the  high-frequency  region  by  only  a  few  terms  of  an 
expansion.  It  has  been  submitted  to  Radio  Science. 
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is  an  unfinished  draft  paper  on  eclipse  shadow  bands,  by  J. 
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Solution  for  the  Fourth  liommt  of  Wa^es 
Propagating  in  Random  Media 
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R.G.  Frehlich,  and  Frank  S.  Henyey 


Center  For  Studies  of  Nonlinear  Dynamics^^ 
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Abstract 

A  series  expression  is  developed  for  the  fourth  moment  of  a  beamed  field 
incident  on  a  random  phase  screen  or  an  extended  medium.  The  series  has  a 
symmetry  that  allows  its  first  few  terms  to  generate  useful  approximations  at 
both  low  and  high  spatial  frequency.  The  parabolic  wave  equation,  the  Markov 
approximation,  and  Gaussian  refractive  index  fluctuations  are  assumed.  The 
result  for  the  phase  screen  is  obtained  by  Green’s-function  techniques.  The 
extended-medium  result  is  derived  in  an  analogous  manner  using  path  integral 
methods.  The  same  results  are  also  derived  by  moment-equation  methods.  The 
behavior  of  the  leading  terms  is  compared  to  previous  results  for  plane-wave  and 
point-source  geometries. 

PACS  numbers:  03.40.Kf.  05.40. 'fj,  42.20.Cc.  92.10. Vz 
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Many  years  of  research  have  been  devoted  to  the  study  of  wave  propaf  ation  in  random 
media  (WPRM).  The  first  comprehensive  review  of  the  field  was  by  Tatarskii  [1971],  followed  by 
Prokhorov  el  al  [1975],  Ishimaru  [1978],  and  Fante  [1975,1960].  The  propagation  of  radio 
waves  through  the  ionosphere  is  reviewed  by  Yeh  and  liu  [1982].  The  phenomenon  of  interstel¬ 
lar  scintillation  is  reviewed  by  Rickett  [1977]  and  Rickett  et  al  [1984].  Sound  propagation 
through  the  ocean  and  path  integral  techniques  are  discussed  by  Flatte  et  al  [1979]  and  Flatte 
[1983]. 

We  consider  waves  propagating  from  an  arbitrary  source  distribution  in  a  random 
medium.  We  assume  the  statistics  of  the  medium  are  locally  homogeneous,  and  sre  make  the 
Markov  approximation:  i.e.  the  field  fluctuations  induced  within  a  correlation  length  along  the 
propagation  direction  are  weak.  For  a  more  complete  discussion  see  Codona  et  al,  [1985]. 
The  wave  propagation  is  characterized  by  narrow  angular  scattering  due  to  the  small  random 
fluctuations  in  refractive  index.  It  is  then  convenient  to  write  the  complex  monochromatic 
scalar  field  as  £(i,z)e^  where  z  is  the  propagation  direction,  lis  the  transverse  coordinate 
and  k  is  the  wavenumber  of  the  wave  with  no  refractive  index  fluctuations. 

The  random  nature  of  the  fields  is  conveniently  described  by  statistical  moments 
evaluated  in  the  transverse  plane  located  at  distance  R.  Ensemble  averages  of  randmn  vari¬ 
ables  are  denoted  by  <>.  The  first  moment 

=<£•(!,/?)>  (1) 

or  average  of  the  field  and  the  second  moment 

ra(l„fB,.R)  =  (le,R)>  (2) 

<»*  mutual  coherence  function  are  well  understood  [Tatarskii,  1971].  However,  there  are  few 
analytic  resiilts  for  the  fourth  moment 


r4(l„le,la,t4./?)  =  {tt,R)E{U.^)^  {K^)> 


(3) 
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Previous  thcoreticol  work  concentrated  on  plane^wave  and  pointfource  geometry.  We 
present  tbree  main  results  for  arbitrary  source  distribution. 

A  series  expression  for  the  fouKh  moment  is  derived  as  an  expansion  of  the  Creen's  func* 
Uon  for  the  fourth  moment,  thus  avoiding  the  difficulties  associated  with  the  source  distribu¬ 
tion.  For  the  thin-screen  problem,  the  expansion  quantity  is  a  combination  of  phase  structure 
functions.  For  the  extended  random  media,  the  expansion  quantity  is  an  analogous  combina¬ 
tion  of  phase  structure  function  densities.  The  Green's  function  is  expressed  as  a  multiple 
path  integral.  The  resulting  series  of  path  integrals  is  evaluated  with  a  useful  identity. 

Our  second  result  is  the  generatiixi  of  two  series  for  the  intensity  correlation  or  intensity 
spectrum.  Ihe  fourth  moment  has  the  obvious  symmetries  that  it  is 

unchanged  by  interchanging  ti  and  ^  or  by  interchanging  ^  and  t|.  Each  term  of  the  series 
expansion  does  not  share  the  symmetry  of  the  entire  expression.  Thus  two  separate  series 
are  obtained  by  invoking  symmetry.  In  principle,  either  series  could  be  summed  to  give  We 
demonstrate,  however,  that  it  is  better  to  cixisider  both  series  in  order  to  describe  the  fourth 
moment  with  the  fewest  number  of  terms.  This  assertion  is  demonstrated  for  the  second 
moment  of  intensity  or  intensity  correlation.  C(tt,^,^).  which  is  a  special  case  of  the  fourth 
moment.  i.e. 

)/(lij./?)>  =  r4(l,.t,.le.i^  =  r4(l,.lalbti)  (4) 

Note  that  the  symmetry  of  the  fourth  moment  has  been  explicitly  indicated.  A  clear  presenta¬ 
tion  of  the  behavior  of  the  intensity  correlation  series  obtained  hrom  the  fourth  moment 
expansion  requires  the  introduction  of  a  spatial  spectrum  of  intensity  fluctuatitms  for  a  spa¬ 
tially  nonstationary  random  process.  We  adc^t  the  definition 

♦(A.4^)  =  dfi  (5) 

where 

A=7(*i  +  ^)  ^  =  (6) 

(Note  the  free  format  of  the  argument  list  of  functions).  The  spectrum  has  the  property 
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(7) 


It  should  be  noted  that  the  spatial  spectrum  may  depend  on  the  centroid  ^ 

Since  there  are  two  series  for  the  intensity  correlation  there  are  also  two  series  for  the 
intensity  spectrum.  The  leading  terms  of  one  series  for  4(/l.4^)  describe  the  small  4 
behavior  while  the  other  series  is  valid  at  high  4  The  rate  of  convergence  of  each  series  pro> 
vides  a  criterion  for  merging  the  two  results  to  produce  a  complete  expression  for  the  inten* 
sity  spectrum.  In  general,  an  analogous  treatment  of  the  intensity  correlation  series  is  not 
possible  since  the  leading  terms  of  both  series  do  not  converge  to  the  variance  as  the  spatial 
separation  approaches  zera 

Our  third  result  is  the  demonstration  of  the  equivalence  of  path  integral  and  moment- 
equation  methods.  Early  theoretical  work  on  WPRlf  concentrated  on  geometrical  optics  and 
the  method  of  small  perturbations  [Barabanenkov.  1971;  Tatarskii.  1971].  These  tsro 
approaches  were  limited  to  weak  scattering  conditiims.  This  restriction  was  removed  with  the 
introduction  of  differential  equations  for  the  moments  of  the  field  [Prokhorov,  1975).  Func¬ 
tional  techniques  of  high  energy  physics  (path  integrals  and  operator  methods)  provided 
another  point  of  view  to  WPRM  [Klyatsldn,  1973;  Dashen,  1979].  The  moment  equation  method 
and  functional  techniques  are  equivalent  [Codona  et  al.  1965]  and  must  generate  identical 
results  when  expansions  are  performed  in  the  same  quantity.  This  equivalence  is  demon¬ 
strated  by  deriving  the  same  fourth  moment  series  expression  using  moment-equation 
methods. 

The  thin-screen  case  is  considered  in  section  2.  The  second  and  fourth  moment  are 
analyzed  with  Green's  function  techniques  and  the  behavior  the  intensity  cmrelation  is 
investigated.  The  same  analysis  for  the  extended  medium  case  is  presented  in  section  3.  Here 
we  use  the  path  integral  representation  of  the  Green’s  function.  Identical  results  for  the 
fourth  moment  are  derived  with  moment  equation  methods  in  section  4.  The  main  results  of 
the  paper  are  summarized  in  section  5. 


^  GREDf ‘S  rUNCnON  itfTIKMCH  TO  THE  THIN  SCROM  PIDBIZII 


2.1  iDtroducUeo 

One  of  the  Arst  WPRM  probiems  considered  was  the  propagation  of  plane  waves  through  a 
random  phase  screen  [Mercier,  1962:  Salpeter.  1967;  Bramley,  1967;  Cocbelashvily  and 
Shishov.  1971.  1972.  1975:  Rumsey.  1975:  Rino,  1979a.  b:  Uscinski  and  Macaskill.  1963a.  b]. 
The  propagation  of  radio  waves  through  the  ionosphere  and  the  solar  wind  are  two  applica* 
tions  of  this  model.  The  theory  of  scintillation  from  a  point  source  viewed  through  a  random 
phase  screen  has  been  investigated  by  Lee  [1977].  The  case  of  a  Gaussian  beam  focussed  on 
the  observation  plane  has  been  considered  by  Cocbelashvily  [1974].  Previous  work  concen¬ 
trated  on  plane-wave  and  pc^t-source  geometries.  We  anal3rze  the  more  general  problem  of 
an  arbitrary  beam  incident  on  a  phase  screen  using  Green’s  function  methods.  The  following 
analysis  is  presented  in  a  fashion  that  permits  a  clear  extension  to  the  more  complex  problem 
of  wave  propagation  in  extended  random  media.  We  review  Green’s  function  methods  with  a 
discussion  of  the  second  moment  A  series  expression  for  the  fourth  moment  is  presented  as 
an  expansion  of  the  Green’s  function  for  the  fourth  moment.  The  behavior  of  the  resulting 
series  for  the  intensity  cwrelation  is  then  discussed. 

Consider  the  scalar  wave  field,  ^(f.s).  incident  on  a  thin  random  phase  screen  sitiuted 
at  the  plane  z=0.  The  field,  ).  emerging  from  this  screen  is  given  by 

^(t'.0  +  )  =  £’(l'.0)c‘*<»>  (8) 

where  £(i!.0)  is  the  field  just  before  the  interface  of  the  screen  and  the  phase  fluctuations  are 

e(i')=*/n(r.*)(fa  (9) 

0 

where  n(f^,e)  is  the  random  fluctuations  in  refractive  index.  Assume  8(t*)  is  a  zero  mean 
Gaussian  random  variable  with  homogeneous  statistics  and  correlation  function 


c.(i)  =  ♦!)>«/ 


(10) 
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whert  is  the  spectrum  of  phase  fluctuations.  The  structure  function  of  phase  fluciua* 
lions.  Pe(|).  is  defined  by 


D^h  =  <[0(t) -e(f  ♦  4)1*>  =  2(Ce(0) - C,(|)] ^zf[\  -cos(A.4)J*e(«  << A  (ID 


For  narrow  angular  scattering  the  scalar  field  sattsfles  the  parabolic  wave  equation 


2ifc|^+  V*£=3 


The  solution  of  the  field  at  a  distance  R  from  the  screen  is 


(12) 


E(%R)  =  f  E{r,0)C{%t.R)  dt 


(13) 


where  the  Green's  function  is 


C(l;t'./?)  =  e‘«{»>C'(l;l'./?) 


(14) 


with  the  free  space  Green’s  function 


(15) 


2.2  Green's-Ftinction  Approach  to  the  Second  Moment 

Since  the  random  fields  have  a  Green's  function  solutioa  the  monents  of  the  field  also 
have  a  Green’s  function  representation.  Consider  the  second  moment  of  the  field  or  mutual 
coherence  function,  r^ii,^.R),  given  by 

m  m 

T^ti.U.R)  ~  <E(KR)^  (UR)>  =  J  {it:tt.R)>  df,  dr,  (16) 

The  Green's  function  for  the  second  moment  is  identified  as  the  expression  inside  the  <>  .  Le. 


(h(U,ti.R)  *  a d,:!-*./? XexpCiCAd*,) -  A(r,))l> 


(17) 
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e( (*,;», =  Jj4^„p{^((t,.r,)*-(k-r,)*))  (js) 

is  the  free  space  Green's  function  for  the  second  moment.  Here,  and  ti  denote  the  set  of  t 
and  t  coordinates  respectively.  The  expectation  over  the  random  phase  is  performed  with  the 
identity 


sHttch  is  valid  for  zero  mean  Gaussian  random  variables.  The  Green’s  function  for  the  second 
moment  is  then 


Chi%:ti.R)  =  Gi  (*i:»*./?)exp[  -  fa)) 


and  the  second  moment  is  given  by 


2.3  Green’e-fnnctioa  Approach  to  the  Vburth  Moment 

Previous  work  on  the  fourth  moment  has  concentrated  on  plane-wave  conditions 
[Zavorotnyi.  19‘79a].  We  address  the  full  fourth  moment  with  an  arbitrary  source  distribution. 
FoUomng  the  previous  analysis  of  the  second  moment,  the  Green's  function  for  the  fourth 
moment  is  given  by 


-  tz)  +  DziU  -  ti)  +  Dziti  - 1*4)  -  ti)  -  D^tz  -  ^4)  -  ^^ed*,  -  ti)]] 


where  G{{ti;ti.R)\s  the  free  space  Green's  function  for  the  fourth  moment.  It  is  convenient 
to  apply  the  unitary  coordinate  transformation  [Rumsey.  1975] 


2asi|4^ 


2S,sA4|.f.a 

2t=l,-I,-I,^l«  21e  =  d-^-f^a 

2a  =  t,-^<t-ta-t«  2t«  =  d-^^t-a  (23) 

The  set  {A'.^.J’.t') ,  will  denote  the  same  transformation  on  the  coordinate  set  (1‘|.1’2.1^9,1'4). 
The  Green's  function  for  the  fourth  moment  then  becomes 


(24) 


exp 


-jWr  ♦  a)  ♦  ^•(r  -  ♦  *•)  ♦  ^#0*  -a*)  -  ♦  r)  -  -r  )l 


where 

C{(Ii:t'i./?)=  |~((i,-i',)*-(ie-ira)*  +  (t8-f3)*-(t«-l’4)«) 


(2n/?)* 


[(a-a').(a-a')+o-a').(f-r)] 


(25) 


is  the  free  space  Green's  function  for  the  fourth  moment  This  expression  is  intractable,  both 
analytically  and  numerically.  Marians  [1975],  numerically  calculated  the  intensity  spectrum 
of  plane-waves  incident  on  a  two-<limensional  phase  screen. 

The  plane-wave  case  was  considered  by  Zavorotn)ri  et  al  [1977].  They  noted  that  a  combi¬ 
nation  of  structure  functions  was  small  in  the  important  regions  of  integration.  A  series 
expression  was  then  obtained  as  a  Taylor  series  expansion.  For  the  general  case,  we  identify 
that  same  expansion  quantity  as 


« =  -  *  3') + x>s(^  -  3')  -  + r )  -  -r  )l 


(26) 


ss  2  y##(d)[cos(3**f)  -co8(f' •f)]e‘^**  dt 
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W«  will  invesiif  ate  the  behavior  of  9  for  a  structure  function  that  is  power  law  above  tbe 
inner  scale  (e-  The  exponent  in  (24)  is  large  unless  two  of  the  first  four  structure  functions  are 
small.  (The  other  two  can  nearly  cancel  the  last  two  structure  functions).  The  only  way  tbe 
exponent  can  be  small .  while  allowing  the  cancellations,  is  for 


i'  =  0(so)  and  =  0(so)  T  -  0(so) 


(27) 


where  Sg  is  the  field  correlation  distance  defined  by  D^s^)  - 1.  The  remaining  variable  (^  or 
^)  is  typically  of  order  of  tbe  scattering  disk 


(2B) 


where  Og  is  the  width  of  the  angular  spectrum.  When  the  inner  scale  is  larger  than  W,  there  is 
little  scintillation.  Therefore,  we  consider  the  case  Sg«W.  lo«W,  ^  =  0(sg).and  ^  =  0(W).  A 
Taylor  series  expansion  about  3*  reduces  Q  to 


^=7(3'*+r®)v0*) 


(29) 


For  a  power-law  structure  function 


(30) 


Q  becomes 


, ,  A), = 0(^-» 


(31) 


On  the  other  hand,  the  other  terms  in  the  exponent  are  0{D^s^)  =  0(1),  which  was  the  condi¬ 
tion  which  caused  f  and  3*  to  be  order  Sg.  Thus,  if  Sg<  W,  a  Taylor  series  in  0  is  appropriate, 
but  a  Taylor  series  in  the  entire  exponent  requires  many  more  terms. 

The  other  possibility,  ^  -  0(sg),  and  f  *  -  0(  W)  requires  an  expansion  in  another  variable, 
O',  obtained  from  Q  by  interchanging  with  f .  This  alternate  expansion  is  the  fundamental 
reason  that  two  different  series  are  required. 
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We  now  return  to  the  Teylor  series  expansion  in  Q,  with  the  result 

s  t  )  =  (H  ♦  1)  ♦  W  - 1)]]  (32) 

naO  * 

I-*-  £  fr/  <«i  f  «<^e*P  '£*,•<**  rt*#(*;)[cos(3'.«,)-cos(y.£,)J 

nai  -•  I 

This  series  should  converge  quickly  when  the  quantity  Q  is  small  over  the  important  region  of 
integration.  Note  that  the  symmetry  of  the  full  fourth  moment  expression  (24)  does  not  hold 
for  each  term  of  the  expansion.  The  equivalent  moment-equation  derivation  is  presented  in 
section  4. 1. 

The  first  term  of  (32)  reduces  to 

=  G2(t,.te.l',.ira/f)Ct(l9.t*.l's.l'4./?)  (33) 

Thus,  for  spatially  coherent  sources,  the  flrst  term  of  the  fourth  moment  is  a  product  of  two 
second  moments,  i.e. 

r^{ti.uuKfi)^W^,%i.R)r2{i^UR)  (34) 

The  next  term  of  the  Green's  function  expansion  is 
G^{ti.t,.R)  =  2G{  (ti:t'„^)exp[ +  «')  +  X?e(r  -«')]]  (35) 

/♦s(£)[cos(J'.£)  -cos(r  •*)]e‘^’*  di 

These  two  terms  contain  the  useful  first-order  description  of  the  foxirth  moment.  The  rate  of 
convergence  of  the  series  is  determined  by  the  higher  order  terms.  For  the  plane-wave  case, 
the  fourth-moment  expressions  generated  from  (34)  and  (35)  are  identical  to  the  asymptotic 
results  oi  Zavorotnyi  [  19*79a  Eq  (7)]. 


Z.A  InUntf ij  CorralaUoo 

There  ere  few  measurements  of  the  full  fourth  moment  of  WPRM  [Curvich  et  al.  1978, 
1979a].  However,  the  intensity  correlation,  a  special  case  of  the  fourth  moment,  is  commonly 
observed.  We  will  now  demonstrate  that  our  fourth-moment  series  (32)  generates  two  different 
expressions  for  the  intensity  correlation  (4),  one  valid  at  low  spatial  frequency,  the  other  valid 
at  high  spatial  frequency. 

The  low  frequency  version  of  the  intensity  correlation  is  obtained  from  the  fourth 
moment  (3)  by  setting  t|  =  ^  and  ^  =  t#  or  by  setting  2  =  0  and  f  =  0.  The  n=0  term  produces 

cVidM  =  ^—^7  r4(d*,^r.i'.o)exp(  -i|-[(d  (36) 

exp[  -  ^[Z)«(r  +  «■)  +  W  -2')]]  dd-  dr  d2* 

2 

For  spatially  coherent  sources 

C\f  =  </(^-|^./?)></(  ^.y?)>  =  </(l,.*  )></(l3.e)>  (37) 

or 


*y(A  Afi)  =  </0i+  f-.i?)></(A-  f  .-ff)>  dfi  (38) 

which  normally  describes  the  smallest  spectral  scale  of  the  process.  Indeed,  for  plane-wave 
conditions  i^(A4-^)  =  </>‘2(4)  since  the  average  intensity,  </>,  is  a  constant 


The  corresponding  terms  for  n=  1  may  be  written  as 


cY(a.M)  =  ^^^/r4(a*>.r.i'.o)e«p 


(») 

exp  -  •^[P»(r  ♦i*)  ♦  i^#(r  -*‘)1  ♦#(«)[coe(l‘.a)  -coe(r**)le‘^**  da*  <£$•  df  dl  di 

and 

♦i/(a  A«)=  (5§!^/r,(a\>’.r  =  -a|-.J-.o)eip[-it-(a-»’).J-le‘»-«-»^  (w) 

exp  -•^(/)e{4^  +  i')  +  ^e(4^“i')J*e(*)lco»(^'**)“eoe(T**)l  aa'dj’  di' di 

For  the  case  of  plane  waves.  (40)  reduces  to  the  familiar  result 

♦r(A4/?)  =  <*a(4)exp[  -/?e(4|-)]8in*Ig*^)  (41) 

which  also  describes  the  low  q  behavior  of  The  Bom,  or  weak-scattering  approxima¬ 

tion  is  obtained  by  ignoring  the  exponential  term.  Similarly,  the  Bom  approximation  for  the 
general  case  is  obtained  from  (40)  by  ignoring  the  last  exponential  term. 

The  high  frequency  version  of  the  intensity  correlation  is  obtained  from  the  fourth 
moment  (3)  by  setting  t|  =  i*  and  Sj;  =  lis  or  by  setting  3  s  0  and  f  =  0.  The  n=0  term  produces 

cy(a.f.y?)  =  ^^|^/r4(a*.3’.r.«'.o>exp  -i|-[(a-a').3'+(f-t).^']  (42) 

exp  -  j[D^f  4  J')  +  Dnif-i’)]  dA’  df  df  dl’  dt 
For  spatially  coherent  sources 

cy  (t,.i8,i?) = r,(i,.«^/?)r;  (t,.ie./f )  (43) 


This  expression  is  the  high  frequency  approximaUon,  i.e.  the  intensity  correlation  is  the 
square  of  the  mutual  coherence  funcUoa  and  may  also  be  derived  by  assuming  that  the 


55 


complex  electric  fleldx  ere  xero  mean  Ceussien  random  vartablei. 
The  nsl  term  reducea  to 


(44) 


exp 


«t«)[cos(3‘.<)-cos(t.2)]e^‘*  dd' df  dl  dt 


and 


Wa.4i?)=  -4f .t.2‘.0)exp(-t|-(a-a')-J')e‘“^**^’''  (45) 


exp 


^c(lt)[cOs(^*l()~COs(f'*lK)] 


Thus  we  have  shown  how  the  fourth-moment  series  generates  two  expressions  for  the  intensity 
correlation:  one  (S  =  O.f  -  0)  useful  at  low  spatial  frequencies,  and  the  other  (t  -  0,^  ~  0)  useful 
at  high  spatial  frequencies.  The  region  of  validity  depends  on  the  statistics  of  the  phase 
fluctuations  and  the  initial  source  distribution.  The  case  of  plane  waves  incident  on  a  random 
phase  screen  with  a  power-law  spectrum  has  been  investigated  by  Gochelashvily  and  Shishov 
[1975].  Their  calculations  of  the  first  few  terms  of  the  intensity  spectrum  [Figure  l]  imply  the 
two  series  converge  quickly  when  the  quantity  Q  is  small  over  the  important  region  of  integra¬ 
tion.  The  rate  of  convergence  is  difficult  to  determine  a  priori;  the  contribution  to  the  inten¬ 
sity  spectrum  from  the  higher  terms  is  the  best  indication  of  convergence.  The  intensity 
correlation  series  is  more  difficult  to  interpret  since  the  ernx*s  of  the  expansion  accumulate 
in  the  region  of  small  spacing  but  in  the  spectral  domain,  these  errors  appear  in  the  central 
regions  of  the  spectrum  [cf  Figure  l]. 
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a  PATH  INITGIULnCHNlQUES  FDR  DCTBIDEDRMfm 

ai  Introduction 

We  now  consider  the  more  complex  problem  of  wave  propagation  in  a  random  media  that 
is  locally  homogeneous  with  statistics  that  vary  slowly  in  the  direction  of  propagatioa  Laser 
propagation  in  the  atmosphere,  radio  propagation  through  the  interstellar  medium,  and 
sound  propagation  through  the  ocean  with  no  deterministic  background  are  common  exam* 
pies  of  this  phenomena.  The  application  of  path  integral  techniques  to  problems  of  wave  pro¬ 
pagation  in  random  media  was  introduced  by  Klyatsidn  and  Tatarskii  [1970].  Zavorotnyi  et  al 
[1977]  and  Dasben  [1979].  Path  integral  methods  have  been  successfully  applied  to  the 
difficult  problem  of  WPRH  for  anisotropic,  inhomogeneous  medium  with  deterministic  back¬ 
ground  of  refractive  index  [Flatte  et  al.  1979].  A  functional  operator  form  of  the  path  integral 
was  used  by  Tatarskii  and  Zavorotnyi  [1980].  to  extend  thin-screen  analysis  to  the  problem  of 
wave  propagation  in  extended  random  media  for  the  plane-wave  case.  We  use  the  path  integral 
representation  of  the  Green's  function  to  illustrate  this  connection  because  one  obtains  a 
clear  presentation  of  the  role  of  the  source  distribution.  Note  that  the  operator  formalism  is 
eqmvalent  to  the  path  integral  method  and  produces  the  same  results. 

The  path  integral  technique  is  introduced  by  a  review  of  the  second-moment  derivation. 
The  reduction  of  path  integrals  to  familiar  Reimann  integrals  is  perfixmed  by  a  useful  identity 
[cf  (59)].  Using  this  identity,  we  present  a  series  expression  for  the  fourth  mcMnent  that  is 
analogous  to  the  thin-screen  results  of  section  2.2.  The  behavior  the  resulting  correlation 
series  is  then  discussed. 

For  narrow  angular  scattering,  the  scalar  field  satisfies  the  parabolic  wave  equation 

2il:~  +  V*£'  +  2fc«n(l;*)£'  =  0  (48) 
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where  n{t,z)  denotes  the  refractive  index  fluctuations.  We  define  the  correlation  of  refractive 
index  fluctuations.  as 


9»(at.t)s<n(0.«)n(M  ♦<)> 


(47) 


Ihc  path  integral  formulation  for  the  Green's  function  was  developed  by  Feynman  [1948],  and 
may  be  written  as 


C(r.l/?)=  /  Pf(s)exp 


i|-/(f(a)]*d*  -at/n(f(s).s]ds 


2  “-(j  -  g 


(48) 


where  Dt{*)  denotes  the  infinite  dimensional  integration  over  all  possible  paths,  F(*)>  con¬ 
necting  the  points  (f.O)  and  (t,R)  and  K* )  * 

If  there  are  no  refractive  index  fluctuations,  the  Green's  function  becomes  the  free  space 

4 

« 

Green’s  function. 


Gf{%t,R)==  f  Dfi*)exp 


« 


ZiriR 


exp 


2/?'* 


(49) 


3.2  Second  Moment  Iqr  Path  Integral  Techniques 

The  Green's  function  for  the  moments  of  the  field  ve  easily  expressed  in  terms  of  the 
path  integral.  The  Green's  function  for  the  first  and  second  moment  were  derived  by  Dashen 
[1979].  These  results  were  obtained  under  the  Mar'tov  approxunation  [Zavorotnyi,  1978].  We 
review  the  derivation  for  the  Green's  function  for  the  second  moment,  which  is  given  by 


G^tx,U±x.Rz,R)  =  <G{tx-,tx.R)d  {Utz.R)> 


(50) 


s  //  Dfxiz )  DM*  )exp 


't/  ((*'(• 


<exp 


-ifcjr[n(f,(s).s  )  -n(fs(s  ).s  )]<le 


Applying  (19)  and  the  Markov  approximation  results  in 
fc* 


«./?)  = 


(2ni?)* 


exp 


fS  Dtit)DUt) 


1 


(51) 
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•nd  the  phase  structure  function  density  is  given  by 


d(|.f ) .  2k*  f  .s )  -  fiU  (I.* .» )ldf 


=  4n4*/[l  -cos(4.^]**(4g,  =  0.s)  d^ 


where 


♦»(49«.^)=  ^^^//ft{t*.^)exp[-i(i.4+sg,)]  dldt 


is  the  refractive  index  spectrum.  Change  path  variables  to 


A(* )  =  7[fi(* )  +  fa(* )]  -  j(l,  ♦  *s)  -I"  -  ji^x  +  i'sXl  - 1-) 


where 


f(* ) = (»I  “  te)  ~  *  (Ir,  -ir*)(i  -  ^) 


Then  A<0)  =  A(/?)  =  Ko)  =  ?(/?)  =  0  and 


C,(l*;r*.y?)  =  exp|  ^[(r ,  -l,)* - f  Dp{z )  pfi* ) 


z)dz  -^J'd[^{z)  -t-7(s).s]iic 


Integrating  the  first  path  term  in  the  exponential  by  parts  and  substituting  the  free  space- 
Green’s  function  results  in 


(88) 


exp 


0  *  0 


This  path  integral  is  evaluated  with  the  identity 


ff  Dp{t)  /)3{*)/'(3(x)lexp 


(2n^)* 


niC*))  (59) 


where  |(x)  is  the  solution  of  i(s)  =  P(e)  that  satisfies  the  appropriate  boundary  conditions. 
Mentity  (58)  reduces  (56)  to 


<.f?)  =  O' (l,:l'i./?)C^  (lfe:ir8.i?)exp 


)••!*** 

•  fi 


(80) 


and  i(s)  =  0  is  the  solution  of  {(s )  -  0  .  that  satisfies  the  required  boundary  conditions.  Using 
thiat  Green's  function,  the  expression  for  the  second  moment  becomes 


r,(i„ie.f?)= 77r2(r,.i'*.0)G,(i,.i^i',.r8.f?)  dr,  dr* 


(61) 


For  plane*wave  input 


r8(l,.fiB.f?)  =  </>e3q» 


(82] 


and  a  point  source  gives 


r2(ti.le./f)  =  </(/?)>exp 


H-w  -*Ji 


(o; 
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3.3  fourth  Momooi  bf  Pith  bUgnl  Tochniqiiot 

Wt  derive  •  series  expression  for  the  fourth  moment  in  e  menner  enelogous  to  the  thin* 
screen  derivation.  The  Green’s  function  for  the  fourth  moment  is  given  by  the  multiple  path 
integral 

<G{tx.ti.F)(f  {i9,tz,R)G{it;Ky.R)d  it^.t^,R)>-  f  Df^  Dtt  Df^  DU 

H  \  M 

«*P  - ^8  +  ^3  - ^41***  <e*p(  ) ~n(r,.* )  +  n(f3.* )  -n(p4.s  )]dx  ]>  (64) 

,0  0 

n 

If  r[n(f|,s )  »n(f2.s  )  4- n(f3,s )  -n(f4,s )]<lx  is  a  zero  mean  Gaussian  random  variable  (19)  and 
*0 

the  Markov  approximation  is  valid.  the  last  term  of  (64)  becomes 
B 

exp[  -  7  /(d(f ,  -  fa.* )  +  ‘f(f3  -f4.* )  MU  -  f4.* )  ♦  MU  -Ut)-  Mh^K* )  -  «*(f  i  -  fs.*  )]<**  ]• 
Change  path  variables  to 

2a(* )  =  f  ,(s  )  ♦  f2(s )  +  f8(*  )  ♦  f4(*  )  2f  i(*  )  =  «(*)  +  )  +  f(*  )  +  ^(*  ) 

2|(t)  =  f,(s)  +  fj(s)-fij(s)-f4(*)  2f*(s)sa(s)  +  |(z)-f(z  )-<(*) 

2T(* )  =  f  »(* )  -  f|(* )  -  f8(* )  +  U(*  )  2f3(*  )  =  «(*)-  ?(* )  - 1(* )  +  H*  ) 

2i(s)  =  f,;z)-^(t)  +  f5(s)-f4(*)  2f4(s)  =  a(z)-^(s)+f(*)-J(z)  (65) 

Then 

Gt{ti:ri.R)  -fDdD0  Di  exp 
B 

exp  /(d(^  +  3.s)4d(^-J,s)  +  d(4  +  J.*)  +  d(^-8,s)-d(^+f.s)-d(^-f,s)](l*  (66) 

I  *  *0 

Change  variables  again  to  paths  about  the  straight  lines  connecting  the  points  (f^.O)  to  (l^.^) 


i.e. 


where 


a,{t  )=«(•) -!(«.#) 


J,(*)  =  3(*  )-!(<,*) 


Ka..) = a(o)(i  -  ♦  4(* )  ^ ♦  a  I"  (es) 

is  the  straight  line  path  from  (d'.O)  to  (a.f?).  In  order  to  simplify  the  large  expressions, 
[a*.#'.t'.3’]  will  replace  [a(0).|(0).9(0).3(0)]  and  [a./l.f.3]  wUl  replace  [a(/?).^(/?).f(/?).3(/?)] 
Then 


C4(V.fi./?)=  DAi  Dh  /?3,exp 


I  . 


expj""”  /[a(Ti  +  3i  +  i(T.3).s)  +  d(fi  -  3i  +  6(^,3).*)  +  d(^i  +  3|  +  a(4.3),s ) 


+  d(^,  -  3,  +  6(^.3).* )  -d(4, + f ,  +  i0.f).s )  -d(^,  -f ,  +  M.f).*  )]d* 


Ka.4)=3(a.s)+f(|.s) 

6<a.4)=3(a.*)“3(^.*) 


Substituting  (integration  by  parts) 
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•nd  performing  Ihc  inUgroUon  over  the  paths  d|  end  l|  using  (SO)  results  in 

H 

e*P  -  I  '*■  I  ♦  .i).» )  +  <<(4i  +  )  (72) 

(i(3,  ♦  )  -  <<(4i  ♦  f  I  ♦  )  -  <<  (^1  -  ti  ♦  6(4.f  ).*  )]«toj 

Note  that  (72)  is  similar  to  (24)  of  the  thin*screen  derivatioa  The  analogous  expansion  param¬ 
eter  is  identified  as  the  combination  of  phase  structure  function  densities 

Q  -  +  *(^.^).* )  +  <i(fi  +  6(^.3).* ) -d(0i  +  f ,  +  lO.f ).*)  -<f0,  -f ,  +  60.f).*)]de 

»  - 

=  4n4*//*n(49«  =0.*)exp[i^.[4i  +  ?(^.x)]][cos(4.F(3.*))-cos(4.[f,  +  ?(f.*)])J  d^dz  (73) 
0  — 

For  a  point  source,  this  expression  is  small  in  the  contributing  regions  of  path  space  [Dashen. 
1979).  Performing  the  Taylor  series  in  Q  results  in 

n 

C4(t*;r,.J?)=  C((li:ir*./?)//  D$,  fJ^exp  -^/[^(f  ,(*) +  i(r«).») +  «i(T.(*)  +  6(r3).*)l«i* 

1*0 

Lxp[-a:/4i(*)«!Fi(* )<!*]+  £  /<**»/ d4i  fd^ 

I  0  mal  0  0  —  -• 

exp[-i*/3,(*).[f,(*)-2  +i£?(^.*/)*^] 

f|  ♦»  (^  .9.  =  0.x,  )[cos  (4, .  f(J.x,  ))  -  cos  (4, .  (Ti(x, )  +  f(T,x, )))]  (74) 

Applying  identity  (59)  yields 


% 
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s  (H (t,;ir,.i?)[eKp[-J-|[(l(i(t J).f )  ♦  d(6(t.l).«  )]dM 


7d^.7d4  /d*. 


ASl 


(75) 


exp 


jm\  *  0  ■»«!  *  * 


=  Z  G^iK.tuR) 

»«0 


where  /i(e:Zi)  is  the  solution  of 


(76) 


that  is 


h(*:*,)  =  *(^-l)  *<*i 


Ai(*;2,)  =  *,(^-1)  *>*, 


(77) 


The  free  space  Green’s  function  for  the  fourth  moment.  G{  (ti'Si.R).  is  given  by  (25). 
The  first  term  of  (75) 


c.  -f7ieWf.D.*)*<CKf.D.*)l<i 

G^JiU.ti.R)^G({ii-.ti.R)e  *•  =  G8(l,.lfe;l',.l'8./?)C8(f^l,:rs.l'4./?)  (78) 

is  analogous  to  the  thin-screen  result  (33). 

^  The  n=  1  term  is  more  complex. 

X  -  - 

G^{ti:ti,R)  -  4n**C( (li :»*.>?)//♦, (4 9.  =  0.*  ,)[cos(4.f(3.*  i))  -cos(4.[  J-h(si.e,)  +  f(^.* »)])] 

0  ^  * 


e*'-’^^'W[-7|[‘*(f^(*:*»)  +  i(f.i).*)  +  <<(f^(*:*»)  +  6(r3).*)l«<« 


d^ds,  (79) 


The  higher  order  terms  of  the  fourth-moment  Green's  function  are  obtained  from  (7b) 
but  become  more  intractable.  However,  the  first  order  description  of  the  fourth  moment  is 
given  by  the  leading  terms  of  the  series  expansion.  We  now  compare  these  expressions  to  pre¬ 
vious  results  for  the  fourth  moment.  The  point-source  result  is  obtained  from  the  Green's 
function  by  setting  to  zero.  This  reproduces  the  iterated  series  of  Shishov  [1972  Eq  (13)] 
that  was  derived  using  moment-equation  methods  in  a  spherical  coordinate  s)rstem.  Appling 
the  plane-wave  initial  condition  generates  the  series  expression  of  Shishov  [1971  Eq  (20)]. 
Gurvich  et  al  [1979b]  proved  this  series  was  convergent  It  can  be  shown  that  the  point- 
source  series  is  also  convergent  Applying  the  plane-wave  case  to  (79)  with  d  =  0  and  3  =  0. 
reproduces  the  strong  scattering  results  of  Fante  [1975  Eq  (Bl)]  and  Zavorotnyi  [19*795  Eq 

(3)]. 

3u4  Intenaity  Correlation 

Furutsu  [1972]  derived  expressions  for  the  intensity  correlation  from  a  Gaussian  beam 
propagating  in  a  random  medium  with  a  square-law  structure  function.  This  case  describes 
random  wander  of  a  beam  [Wandzura.  1960].  Weak  scattering  results  for  intensity  correlation 
from  a  Gaussian  beam  are  presented  by  Ishimaru  [1969].  We  now  consider  the  intensity 
correlation  following  the  thin  screen  analysis  of  secUon  2.3. 

The  fourth-moment  series  generates  two  expressions  for  the  intensity  correlation;  one  is 
obtained  by  setting  3  =  0  and  f  =  0  (the  low  spatial  frequency  region)  and  the  other  by  setting 
3  =  0  and  3  =  0  (the  high  spatial  frequency  region).  The  resulting  expressions  for  the  a=0  term 
are 


c\fid,lR)=  ^^^^7r4(d’.3'.r.3'.0)exp[-i|-[(a-d').3'+(3-3')-r]]  (eo) 


R 

exp  -7/[d((r  +  3')(l--|-).*)  +  d((r-3')(l--|-).x)]d*  dd  df'dfdl 


65 


For  spatially  coherent  sources,  these  expressions  reduce  to 


Cy  {d$.R)  r  </(  /?)></(  ^^,R)>  =  </(!,.*  )></(V* )>  (82) 


cy  (ti.ii./?) = r,(i,.«e.j?)r,*  (i,.*e./?)  (es) 

which  are  analogous  to  the  thin-screen  results  (37)  and  (43). 

The  corresponding  expressions  for  the  ns  1  term  are  given  by 

cY(M.f?)=“~j7-|/r4(a'.^*.r.2'.o)c  ^  c  ^  *  ♦n(4.*i) 

exp  -^|[d(fM*;*,)  +  (r  +  i‘)(l-|-).*)  +  ‘<(f-M»:*i)  +  (r-3')(l-;|-).*)l«<* 

[cos(A*«'(l-  ^)-cos(^.( |-h(* ,.*,)  + r(l -—)))]  da-  tip  dfil’d^dz,  (04) 


-ii-Ut-ri-rtF-W-t)]  if-ru-^l 


*  ♦»(4*i) 


exp  -^/[d(J-A(*:*,)  +  (r  +  3')(l-|-)  +  f|-.*)  +  ‘<(f^(*:*i)  +  (r-«’)(l-|r)*f^.*)ld* 
jcos[4*i'(l-^l-cos[V(|-M*».*i)  +  r(l-X^  +  tX^l  dd-d^-dtdJ'd^d.,  (05) 
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The  corresponding  expressions  for  the  intensity  spectrum  ere  obteined  from  (44).  The  Bom 
•pproximetion  is  secured  from  (B4)  by  ignoring  the  lest  exponentiel  term. 

The  tow-frequency  series  for  the  intensity  spectrum  hes  been  investigated  for  a  power-law 
spectrum  of  refractive  index  fluctuations  end  plane-wave  [Cochelashvily  end  Shishov,  1971, 
1974]  and  point-source  [Cochelashvily  et  al.  1974]  geometries.  The  qualitative  behavior  of  the 
intensity  spectrum  is  similar  to  the  thin  screen  result  [Figure  l].  In  weak  scattering,  the  first 
two  terms  describe  the  complete  intensity  spectrum.  In  strong  scattering  the  intensity  spec¬ 
trum  is  characterized  by  two  components.  The  low  frequency  region  is  described  by  f 
The  high  frequency  behavior  was  first  determined  by  summing  the  low-frequency  series,  which 
yields  the  high  frequency  approximation  The  derivation  presented  here  generates 

this  same  result  plus  higher  correction  terms  without  performing  a  complicated  summation. 

Uscinski  [1982]  derived  an  approximate  expression  for  the  intensity  spectrum  for  the 
case  of  plane  waves  by  summing  a  perturbation  series  for  the  fourth  moment  equation. 
Macaslcill  [  1983]  derived  the  same  result  using  the  two-scale  expansion  [Beran  et  al.  1982].  In 
our  notation,  their  expression  for  the  intensity  spectrum  is 

(2^^**** 

-^d(l+ J-(f?-z).*)-^d(l-|-(^-e).z)]d*C“*<-»  dt  (86) 

The  low  and  high  spatial  frequency  behavior  is  identical  to  that  obtained  from  (84)  and  (83). 

The  intensity  correlation  from  a  point  source  embedded  in  a  random  medium  with  irre¬ 
gularities  that  are  constant  in  one  transverse  dimension  is  obtained  from  the  point-source 
result  by  the  substitution  4it($  =  ^ii(?s)^(9y  )>  The  leading  order  terms  of  the  intensity  corre¬ 
lation  series  are 

=  </>*exp[  -|d(^^)dz] 


(87) 
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CY(fi./f)-  8nk*</ >*/ 7 ♦,(g .1  |r-^(*  i.*  0) 


/t  g 

exp  -^<f[^Ai(x:t,).f  ]<li  ■‘■*^9*^  ^ 


ip.R)  =  8Tri:*</ >*/ /♦nCfl  .*i)ain*[9  u*  i)l 


ZR  Zk 


K  , 

exp  *<9  <<*1 


Ignoring  the  exponential  term  of  (68)  reproduces  the  weak  scattering  (%tov  approximation) 
result  of  Tur  and  Beran  [1963]. 

Fante  [1983]  investigated  the  effect  of  the  inner  scale  of  turbulence  on  scintillation  for 
the  case  of  strong  scattering  of  plane  waves  by  calculating  the  expression  C*/  (f,R).  The  con¬ 
tribution  of  this  term  to  the  total  intensity  variance  is  appreciable  over  a  wide  range  of 
parameter  space.  TTierefore.  higher  terms  are  required  for  an  accurate  description  of  the 
intensity  correlation. 

4.  HOKEKT  EQUATION  APPROACH 


4.1  Moment  Equation  Method  for  the  Thin  Screen 

We  now  derive  the  fourth-moment  result  (24).  using  the  moment-equation  method.  This 
method  is  based  on  partial  differential  equations  for  the  transverse  moments  that  are  derived 
from  the  parabolic  equation  for  the  random  6elds.  These  differential  equations  may  then  be 
solved  by  transform  methods  [Rumsey,  1975]  and  the  method  of  characteristics  [Kiang  and 
liu.  1962].  There  are,  in  general,  many  transforms  that  will  simplify  the  problem.  We  choose 
one  [Shishov,  1971]  that  permits  an  analogous  derivation  for  the  fourth  moment  of  waves  pro¬ 


pagating  through  extended  random  media. 


The  fourth  moment  setitfles  the  dilTerenUel  equetion 


(») 


with  initial  condition 


e*p[  -  -  *8)  ♦  ^e(*a  “  t«)  ♦  D^ti  -U)*  ^>#(*8  -  *9) "  ^#(*8  “  *•)  “  ^•{*« '  *9)31  (®1) 


Change  variables  to  the  coordinate  system  (23).  Then 


Tliis  equation  is  solved  in  the  Fourier  transform  domain 


which  transforms  (92)  to 


If  = 


Changing  variables  to 


A-kl~tz 


r  =s 


3(r)  = 


*+«r 


*  sr 


9,=ty, 


^=-4.v, |r 


transforms  (94)  to 


^(«.4i.ar)=o 


(92] 


(93 


(94 


(9! 


(9 


Sinct  M  is  independent  of  r,  the  solution  in  the  orifinel  variables  is 


-  t  f  .0  + ) 


m 


The  fourth  moment  is  the  inverse  Fourier  transform  of  this  expression.  After  some  manipuia' 
tion,  the  Green's  function  result.  (24),  is  produced.  The  moment'  equation  method  is  based  on 
the  formulation  of  the  free  space  Green's  function  in  the  Fourier  transform  domain.  This 
transform  technique  will  now  be  applied  to  the  more  difficult  problem  of  wave  propagation  in 
an  extended  random  medium. 


4.2  Momait  Equation  Method  for  Extended  Random  Media 

The  moment-equation  method  was  a  major  contribution  to  the  theory  of  wave  propaga¬ 
tion  in  extended  random  media  [Prokhorov  et  al  1975.  Tatarskii.  1971].  Using  moment- 
equation  methods,  a  series  expression  for  the  fourth  moment  was  derived  for  plane-wave 
[Shishov,  1971]  and  point-source  [Shishov,  1972]  conditions.  For  an  arbitrary  source  distribu¬ 
tion.  we  present  an  analogous  expression  that  is  identical  to  the  path  integral  resxilts  of  sec¬ 
tion  3.3.  thus  demonstrating  the  equivalence  of  the  two  methods. 

The  transverse  moments  of  the  electric  field  satisfy  differential  equations  [de  Wolf,  1967; 
Brown.  1967;  Shishov.  1966;  Dolin.  1966;  Chernov.  1969;  Tatarskii.  1969;  Lee,  1974].  We  fdlow 
the  techniques  of  Shishov  [1971]  to  solve  the  diflerentied  equation  for  the  fourth  moment 


^[Vf-Vf +  71-74*^4+ Kr4  =  o 


(96) 


where 


2  K  =  d(l,  -  te.*  )  +  <<  (*i  -  »«.* )  +  (le  -  la* )  +  d(J8  - 1«.*  )  -  d(l,  -  1b.*  )  -  d  (le  -  t».s ) 


(99) 


In  the  coordinate  system  of  (23).  the  fourth  moment  satisfies 


ar4 


a7--rlV7|+v,.v,]r4+Kr,  =  o 


a 


(100) 
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with 

We  identify  the  last  four  structure  function  densities  in  K ,  i.e. 

25  =  d(^  +  a.t)  +  d(3-3.s)-d(|-ff.s)-d(4-f.i)  (lOe) 

as  the  related  quantity  to  the  path  integral  expansion  variable  Q  of  (73).  The  Fourier 
transform  (93)  converts  (100)  to  the  integral  equation 

J-.V^  +  -^(d(t+a)>d(t-a))i/  =  C(W.8.s)  (103) 

where  C  is  the  convolution  of  II  with  the  transfcxTn  ot  S,  Le. 


(^.s  )[cos(a .  4,)  -cos(f .  4)]lf  (*4“  Aj-I*  .*  )d^  (104) 


The  change  of  variables  (95)  transforms  (103)  to 


^  ♦  jtiW’-)  ♦  Hr).T)  *  (i(f<r)  -Kr).r)Mf  =  C(t.4«r).J(r).r) 


(105) 


lliis  is  an  ordinary  differential  equaticm  in  r.  with  solution 


MM^i.r)  =  II  p0)exp 


-  jf  [<*  (f(*  )  +  H*  ).»)♦«*  (^*  )  -  a(s  ).s  )]dc 


r 

/exp 


"  (f(* )  +  i(*  ).*)  +  «<  (^*  )  -  H*  ).*  )]<te 


c(d4?(*i).<(*i).*i)  *i  (106) 


Changing  variables  back  to  Z,%R  results  in 


Jl/(«.4f.a.i?)  =  ji/(*.4ts.3s.O)exp 


-  •^/[d(^2  ^8.*  ^8.*  )]<** 


♦  /exp 


“ -/[^(fi  ♦  )  ♦  <<(f8  -  Ja*  ))«<« 


G{i.^%X*x)d*x  (107) 
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where 


ii  =  3o(^“*i)  3a  39=3o(/?) 


(108) 


3o(x)  =  3-|-* 


(109) 


Tliis  expression  can  be  written  as  the  integral  equation 


jr .  .  . 

0 


(110) 


where 


R 

ZitA%lO)  =  A/(«.4^3.3s.0)exp  -7/[<i(ta  +  3a.*)  +  d{%  -  3*.*  )]dz 

*  0 


(111) 


.f  *.3.3'./?.* ,) = 4rr***n  (^-  ^  .*  i)3(3'  -  3j)3(f*  -f ,) 


(112) 


exp  -■J-/[rf(l9  +  3a.»)  +  «i(f9-38.*)]d»  tos(3'.(4-4)] -co5(f  .(4-^)]) 
*« 


The  formal  solution  is  written  as  a  s\im  of  terms  which  are  obtained  by  the  iteration 


R  m  m  m 

Jt{*.4f.3./?)*^///K(«.4^.tf'.3.3'.y?.*,)if*-,(*.«.r.i’.*»)  df  dl  dz. 


(113) 


and  Mo- Z.  'Rus  expression  is  a  power  series  in  the  group  of  phase  structure  function  densi¬ 
ties  denoted  by  Q,  which  is  the  same  quantity  as  (73)  of  the  path  integral  expansion;  hence 


the  two  series  ere  equal,  term  by  term.  Unfortunately,  for  higher  order  terms,  this  equality  is 
not  obvious  by  examination  since  the  functional  form  of  the  two  series  is  very  different:  the 
moment  equation  series  is  essentially  a  multiple  convolution,  while  the  path  integral  result  is 
not.  However,  with  careful  algebraic  manipulation,  the  equality  of  the  two  series  can  be  expli¬ 
citly  demonstrated. 

The  moment-equation  method  was  first  applied  to  the  simple  geometries  of  plane  waves 
and  point  sources.  Since  the  moments  are  then  independent  of  the  centroid  coordinate,  the 
differential  equations  simplify.  Early  work  concentrated  on  moment-equation  methods 
because  of  this  simplification.  Moment-equation  teclmiques  have  also  been  used  to  investigate 
the  validity  of  the  Markov  approximatim  [Klyatskin.  1969;  Klyatsldn  and  Tatarskii,  1971]  and 
the  effects  of  non-gaussian  refractive  index  fiuctuations  under  the  Markov  ai^roximation 
[Klyatskin,  1975].  The  Green’s  function  formulation  provides  a  clear  connection  between  the 
thin  screen  and  extended  medium.  The  operate  form  of  the  path  integral  has  also  been  use¬ 
ful  for  evaluating  the  corrections  to  the  Markov  approximation  for  the  higher  intensity 
moments  [Zavorotnyi,  1978]. 

5.SUiaiART 

A  series  expression  has  been  derived  for  the  fourth  moment  of  waves  incident  on  a  phase 
screen  or  propagating  through  extended  random  media.  The  asymmetric  terms  of  the  expan¬ 
sion  generate  two  expressions  for  the  intensity  correlation;  one  that  approximates  the  low  fre¬ 
quency  region  of  the  spatial  spectnim  and  the  other  approfx'iate  fw  the  high  frequencies. 
Ihe  rate  of  convergence  of  the  two  approximations  can  be  used  to  produce  a  complete 
expression  for  the  intensity  spectrum  valid  for  any  initial  source  distributioa  These  results 
can  be  derived  using  moment-equation  techniques  or  functional  methods  (path  integral  or 
operator). 

Many  problems  of  VPRM  may  now  be  attacked  with  renewed  confidence  since  we  have 
derived  expressions  that  are  valid  for  arbitrary  source  distributions.  These  include 
a)  The  intensity  statistics  from  beamed  lasers,  navigational  beacons,  radar,  spacecraft  and 


MltlliU  radio  transmissions,  compact  astronomical  sources  and  other  extended  *avc 
sources. 

b)  The  effects  of  slowly  varying  refractive  index  statistics  along  the  propagation  path.  e.g.  the 
turbulence  profile  of  planetary  atmospheres. 

c)  Comparison  of  thin-screen  and  extended-media  results. 

d)  The  relationship  between  intensity  sUtistics  and  the  form  of  the  refractive  index  spec¬ 
trum. 

e)  Imaging  through  random  media. 

fidmofoledgementa:  This  research  was  funded  by  DABPA. 
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Figure  1.  The  leading  terms  of  the  intensity  spectrum  versus  normalized  spatial  fre> 

i. 

quency.  qRf,  where  Rj-(R/k)*  is  the  Fresnel  scale.  The  curves  are  calculated 
from  expressions  given  by  Gochelashvily  and  Shishov  [  197S]  for  the  case  of  plane 
waves  incident  on  a  random  phase  screen  with  a  Kolmogorov  spectrum  of  phase 
fluctuaUons  and  D%{Rf)~  100.  The  (*)  sign  indicates  that  is  negative  at 

high  frequency. 


IV.  UNEQUAIaFREQUENCY  CORRELATION 


The  following  section  is  a  draft  paper  on  the  two-frequency  intensity  corre¬ 
lation.  This  subject  is  important  for  understanding  the  maximum  bandwidth 
available  in  a  communications  channel.  It  is  also  valuable  in  preparing  for 
understanding  the  connection  between  pulse  propagation  characteristics  and 
the  intervening  random  medium.  It  is  not  directly  related  to  pulse  width 
because  that  depends  on  the  correlation  of  particular  frequency  derivatives  of 
the  fourth  moment  rather  than  the  fourth  moment  itself. 
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Abstract 

The  intensity  cross-spectrum  (spatial  Fourier  transform  of  the  two  frequency 
intensity  correlation)  for  scintillati(Mis  caused  by  a  plane  wave  passing  through  a  ran¬ 
dom  phase  screen  is  considered  By  examining  the  two  frequency  fourth  moment,  we 
show  that  the  "fully  saturated”  approximation  of  modelling  this  spectrum  as  the 
transform  of  the  square  of  the  second  moment  breaks  down  when  the  outer  scale  is 
large  compared  with  the  diameter  of  the  scattering  disk.  A  new  approximation  is 
obtained  which  has  essentially  the  same  limits  of  validity  as  the  "fully  saturated" 
approximation  to  the  monochromatic  fourth  moment 
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1.  Intreduetloo 

The  finite  bandwidth  of  reel  lources  end  receivers  hes  an  important  effect 
on  the  intensity  variation  of  waves  propagating  in  random  media:  scintillations 
will  appear  weaker  than  with  an  ideal  narrow>band  receiver  due  to  intensity 
decorrelation  at  different  wavelengths  [Budden  k  Uscinsld.  1970;  Fante.  1977], 
In  weak  scattering  conditions,  the  decorrelation  bandwidth  is  larger  implying 
that  the  reduction  of  the  intensity  variance  s  minimal  However,  in  strcxig 
scattering  the  effect  of  this  decorrelation  can  be  substantial  and  is  evident  in 
observations  of  "dynamic  spectra".  e.g.  radio  frequency  spectra  as  a  function  of 
time,  which  are  commonly  used  in  analysis  of  pulsars  [Rickett,  1969].  Ibis  effect 
is  particularly  important  when  observing  the  intensity  spread  of  propagating 
pulses  (pulse  jitter)  [Liu  &  Yeh,  1971]. 

Previously,  intensity  decorrelaticms  were  studied  using  a  "saturated  scin- 
tillatiuis"  argument  [Shishov,  1973;  Lee,  1976].  When  the  intensity  scintillations 
have  saturated,  it  is  argued  that  the  real  and  imaginary  parts  of  the  field 
become  zero-mean  Gaussian  random  variables,  implying  that  the  intensity  dis¬ 
tribution  is  Rayleigh.  The  correlation  of  intensity  fiuctxiations  is  then  the  magni¬ 
tude  squared  of  the  second  moment  of  the  fields.  This  is  reasonably  well  under¬ 
stood  for  the  monochromatic  intensity  correlation  (second  moment  of  inten¬ 
sity),  where  it  can  be  shown  that  the  Raideigh  limit  is  reached  in  strong  enough 
scattering  [Tatarskii.  1971;  Flattd.  1983;  Codona  et  al  1985].  When  there  is  a 
finite  bandwidth  present,  the  major  problem  with  the  argument  that  the  two- 
hrequency  intensity  covariance  is  the  magnitude  squared  of  the  two-frequency 
second  moment  (mutual  coherence  function)  is  the  appearance  of  the  mean- 
square  phase  shift  in  the  result  Physically,  this  is  unreasonable  since  the  inten¬ 
sity  of  the  wave  shouldn't  be  affected  by  its  integrated  phase.  This  dependence 
on  the  phase  shift  also  implies  that  in  situations  where  the  important  region  of 
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scattering  (scattering  disk)  is  much  smaller  than  the  outer  scale  of  the  tur* 
bulent  fluctuations,  the  result  from  the  "saturated  scintillations"  argument 
depend  sensitively  on  this  outer  scab  (through  the  phase  shift  factor).  Previ* 
ously,  this  dilemma  had  been  solved  by  dropping  the  phase  factor  entirely.  In 
this  paper  we  analyze  the  question  intensity  decorrelation  of  a  plane  wave 
propagating  through  a  random  phase  screen,  k  series  representation  of  the 
two-frequency  intensity  correlation  hnction  will  be  derived.  The  nature  of  the 
above  phase  shift  factor  will  be  clariled  and  the  conditions  under  which  it  can 
be  ignored  will  be  presented. 

¥«  consider  plane  waves  nonrnlly  incident  on  a  phase  screen  and  freely 
propagating  a  distance  R  beycmd  the  screen  to  an  observation  plane.  The  refrac¬ 
tive  index  fluctuations.  m(^)«  In  the  screen  induce  an  random  phase  change. 
Osk  f  ,  as  the  field  passes  through  the  screen.  Here  s  is  the  direction 

of  the  propagating  wave,  t  is  the  coordinate  transverse  to  this  direction,  and  k 
is  the  wavenumber.  The  phase  change  9  is  assumed  to  be  a  zero-mean  Gaussian 
random  variable  with  locally  homogeaeous  statistics,  and  the  correlation  func¬ 
tion 

C,(l)  =  <9(l')e(l>l)>  =  /♦,(d)e‘«  »  (1) 

where  >3  the  spectrum  of  phaie  fluctuations  and  <>  denotes  ensemble 

averages.  Since  we  are  not  dealing  wib  monochromatic  fields,  it  is  important  to 
remember  that  there  is  wavenumber  lependence  in  the  phase  and  is  included  in 
the  phase  spectmm.  The  wave  propafation  is  characterized  by  narrow  angular 
scattering  in  the  phase  screen  due  to  the  small  fluctuations  of  the  refractive 
index.  Then  the  complex  scalar  wave  field  can  be  expressed  as  E{Lz\k)m^ , 
implying  that,  away  from  the  screen,  be  field  E  satisfies  the  parabolic  equation  [ 

ar  _ 

8s  "2fc  at« 


(2) 


with  the  value  on  emerg inf  from  the  eereen 

eiXO*  Jt )  *  £<1.0;* )  i‘*»>  (3) 

expressed  in  terms  of  the  incident  field  £(1.0:* ).  For  a  plane  wave  the  incident 
field  is  a  constant,  which  we  set  to  one.  The  solution  of  the  parabolic  equation  at 
the  observation  plane  with  the  proper  initial  condition  is  then 

^(t./?;*)=  ^|^/e**(*^exp(^t^*’)* )  d*x'  (4) 

Averages  of  the  random  phase  are  simplified  by  the  use  of  the  identity 

<e«*>  =  e  “  (5) 

valid  for  any  zero-mean  Gaussian  random  variable.  This,  for  example,  gives  the 
average  ot  the  field  as 

<£(*.£;*)>  =  exp -|4*  (6) 

where  *^-Ce(9)  is  the  mean  square  phase  shift. 

In  Section  2  we  introduce  the  general  second  and  fourth  moments  of  the 
field  and  the  two-point,  two-frequency  (intensity)  correlation  function.  In  Section 
3  we  consider  the  "fully  saturated"  approximation  to  this  correlation  function. 
In  Section  4  two  different  series  representations  for  the  intensity  cross¬ 
spectrum  (spatial  Fourier  transform  of  the  intensity  correlation)  are  presented. 
The  first  series  describes  the  low  spatial  frequency  behaviour  of  the  cross¬ 
spectrum  (Le.  it  converges  quickly  in  this  region).  The  other  series  describes 
the  high  spatial  frequency  behaviour.  These  series  are  very  similar  to  those 
fdiicb  appear  for  the  monochromatic  intensity  correlation.  It  is  shown  that  the 
"fully  saturated"  approximaUon  (with  the  mean-square  phase  shift  factor)  is 
only  valid  over  a  negligible  portion  of  the  intensity  cross-spectrum.  The  condi¬ 
tions  under  which  this  phase  shiR  factor  may  be  omitted  are  detailed.  Finally, 


in  Section  5  we  eccumuleie  the  formulee  which  eppronmete  the  two*frequency 
intensity  cross*spectruin.  This  eppnnimetion  is  sufficient  to  analyw  most  prob¬ 
lems  of  intensity  decorrelation  caused  by  strongly  scattering  thin  screens. 


2.  Sbeond  and  Fburth  Moments 


The  random  nature  of  the  fields  is  conveniently  described  by  statistical 
moments  evaluated  in  the  observation  plane.  The  moments  of  particular  con¬ 
cern  for  the  study  of  frequency  decorrelation  are  the  second  and  fourth 
moments.  Ihe  general  secmid  moment  is  given  by 


rz{±i.±z.R:k  M  s  <£(t,.7?;i:,)£*(ta./?;*s)> 


ktkt 

i2itRy‘ 


Lif 


exp 
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where  we  have  used  the  solution  for  the  fields.  (4).  We  have  used  (5)  to  obtain 
exp[-lV2] .  where 


*,-*2 


(ki^kz)/2 


♦«  +  D{tr-tz) 


(8) 


We  have  expressed  the  quantity  Vz  in  terms  of  which  is  the  mean-square 
phase  shift  for  a  wave  with  the  mean  wave  number  k~{ky-^k^/2,  and  D  is  the 
phase  structure  function  associated  with  a  wavenumber  which  is  the  geometric 
mean,  (fciisa)*^*.  Le. 

It  will  be  useful  to  explicitly  display  the  wavenumber  dependence  of  the  phase 
structure  function  by  expressing  it  in  terms  of  a  structure  function  for  the 
integrated  index  of  refraction  fluctuations.  D,^  via  D-kik^D^  where 


I 


■  <(/ A'f>  (10) 

li  is  convsment  to  express  the  second  moment  in  terms  of  mesn  end  difference 
variables  for  the  co*ordinates  in  the  observation  plane 


(11a) 


with  similar  expressions  for  the  (primed)  co-ordinates  on  the  screen.  We  also 


use  combinations  of  the  wavenumbers  given  by 


its*,-*, 


Then  the  second  moment  is  given  by 


(11b) 


kikz  /6k^ai^  r  r 


exp  i^(d-a')*+0-?')*]+i|^a-a')(?-^*)l  dVci*/r 


We  can  immediately  do  the  a’  integration  since  the  phase  structure  function  is 
independent  of  this  variable.  Note  that  this  implies  that  Pa  is  independent  of  A, 
vdiich  follows  from  the  translation  invariance  of  the  problem  for  an  incident 
plane  wave.  Performing  the  a'  integratimi  we  secure  for  the  second  moment 


£  2 


The  general  fourth  moment  is  given  by 
r4(l,.ie.fs.l4.^:*i.*2.it3.*4)=<^(»i.^:*i)^'(»B.^;fcs)^(*i.^:*3)^*(*..^:*4)> 


=//7/rth 


•exp[- 


(-iy*i 


(*f-*,*)*lk 


where  K4  ii  deflned  as 


^4(»i '.*»)=<[/(*  ')-*4M(ti'.«’))de  •]*> 

(15) 

r 

It  is  convenient  to  transform  the  co-ordinates  of  the  observation  plane  to 

a]  1  1  1  1] 

1-1-1  ±z 

^-21-1-1  1  ij 

a]  1-1  1-1)  [i^ 

(similarly  for  the  primed  co-ordinates  of  the  screen),  and  express  combinations 
of  the  wavenumber  variables  as 


£ 

1 

1 

1 

1 

_  1 

1 

1 

-1 

-1 

t  "4 

1 

-1 

-1 

1 
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Afe 

1 

-1 

1 

-1 

In  terms  of  these  variables  the  full  fourth  moment  is 


The  integration  over  the  centroid,  a'  is  immediate  since  the  functions 
depend  only  on  co-ordinate  differences.  The  resulting  expression  is  relevant  for 
calculations  of  quantities  like  pulse  jitter.  The  quantity  we  are  concerned  with 
tor  frequency  decorrelation  is  a  special  case  of  the  fourth  moment  -  the  two- 
point.  two-frequency  intensity  correlation  function 
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Because  of  translation  invariance,  C/  is  a  function  only  of  the  diflerence  variable 
In  terms  of  the  transformed  variables  in  (IS-l?).  we  And  that 
f-o~^-k-0  so  that,  for  the  intensity  correlation.  (IB)  reduces  (after 
integrating  over  the  a'  variable)  to 

■i-'’*  «(!■*  ^r)  (20) 

(?-?')  txp  d’fd*Yd^f 

where  k-{k^-¥k^/2,  6k-ki-k^.  and  is  the  reduced  form  of  (15) 

v'4{4’.f  ■)=*?  (21) 

The  intensity  cross  spectrum.  4/,  is  given  by  the  ^  Fourier  transform  of  (20). 
Before  proceeding  to  evaluate  the  cross-spectrum  we  will  discuss  the  “saturated 
scintillations"  approximation  to  it 

3.  FXiUy  Sahirated  Apimadmatiaa  fix'  Intenaity  Correlations 

If  the  Beld  is  a  zero-mean  Gaussian  random  process  then  the  intensity 
correlation  function  can  be  expressed  as 

In  the  thin  screen  model  where  we  are  considering  an  incident  plane  wave,  the 
mean  intensities  in  the  Arst  term  are  independent  of  position.  The  second  term 
is  only  a  function  of  as  noted  above.  Inserting  Fi  from  (12)  and 

transforming  over  fi  yields  the  saturated  approximation  to  the  intensity  cross¬ 
spectrum 


\26k  R 


[(l,)Mla)*l  6i0,-h^§^i)d%d^0^ 

*1*3 


Changing  integration  variables  to 


and  defining 


we  secure 


There  are  two  problems  with  the  satiu-ated  approximation:  the  mean- 
square  phase  shift  factor  and  the  low  spatial  frequency  behaviour.  In  mder  to 
see  where  the  difficulties  arise,  consider  a  pure  power-law  structure  function 


D{i)  =  \ 


where  sq  is  the  coherence  length  of  the  medium  (  D{sq)-1).  Define  the  size  of 


the  illuminated  spot  (scattering  disk)  as 


R,^^,R 


where  d,  is  the  width  of  the  angular  spectrum,  and  a  measure  of  the  "strength  of 
scattering" 
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The  other  quantity  of  interest  is  the  outer  scale,  Imim’*  of  spectrum  of  phase 
irregularities  where 

/?(•)=♦*  s>l^„  (30) 

In  most  cases  of  practical  interest,  the  outer  scale  is  much  larger  than  the  size 
of  the  scattering  disk.  In  terms  of  R,  the  saturated  approximation  to  the  intn- 
sity  cross'spectrum  is  written  as 

+  exp-(^«** 

exp  (31) 

The  cross-spectrum  is  negligible  for  wavenumbers  greater  than  the  inverse  of 
the  coherence  scale  so  that  the  combination  soc  appearing  in  (31)  is  smaller 
than  unity  in  the  region  of  interest. 

The  argument  for  the  field  to  be  a  zero-mean  Gaussian  random  process 
begins  by  assuming  that  the  coherence  length,  sq  is  small  compared  to  the  size 
of  the  illuminated  spot.  This  implies  that  the  field  at  the  observation  point  is  the 
sum  of  very  many  independent  contributions  from  the  scattering  disk.  By  the 
centred  limit  theorem  the  fields  appear  to  obey  Gaussian  statistics  and  the 
intensity  follows  a  Rayleigh  distribution.  However,  when  the  outer  scale  is  larger 
than  the  scattering  disk,  there  are  large-scale  phase  fiuctuaticms  which  lead  to 
coherent  contributions  of  the  fields  from  different  regions  of  the  illuminated 
spot,  so  that  the  fields  don't  behave  as  Gaussian  processes.  One  immediate 
consequence  of  this  picture  is  that  there  are  intensity  fluctuations  over  regions 
larger  than  the  scattering  disk.  This  effect  will  modify  the  cross-spectrum  at 
wavenumbers  ksS\/  R,  and  is  described  by  geometrical  optics  at  small  enough 
wavenumber.  How  this  modification  occurs  from  a  fourth-moment  calculation  is 
considered  in  the  next  section. 
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The  randomness  induced  by  the  phase  fluctuations  much  smaller  than  F, 
are  consistent  with  the  fleld  being  a  Gaussian  process.  Thus  the  "saturated 
approximation"  is  roughly  correct  if  we  replace  ^  in  (26)  by  D{R,)  of  (29).  In 
the  next  section  we  show  how  this  procedure  leads  to  a  good  approximation  for 
the  large  x  behaviour  of  the  crosa>apectrum. 


4.  Intenaily  Croaa^iectrum 


The  intensity  cross  spectrum  is  given  by  the  Fourier  transform  of  (20), 


which  just  produces  the  delta  function.  R),  where  t  is  the 

iC 


>-f4) 

transform  variable.  With  this  delta  function.  (20)  reduces  to 

*  —  ^  y 


(32) 


with 


2k 


(33) 


=  kfD^{R,Sot{U^)  +  k§D^iR,sot{l-^) 
2k  2k 


6k 


-MV#-  2k 


Unfortunately,  the  final  integral  in  (32)  cannot  be  done  exactly.  Our  approach 
is  to  find  two  separate  approximations  for  high  and  low  spatial  frequency.  Pu^ 
we  consider  high  spatial  frequencies  because  this  is  the  region  where  the  mean- 
square  phase  shift  term  in  (31)  appears. 

In  the  monochromatic  case  at  high  spatial  frequencies  and  strong 
scattering,  the  terms  in  Ki  that  dominate  the  integral  for  the  cross-spectrum 
are  [Gochelashvily  k  Shishov,  1975;  Rumsey,  1975]. 

Separate  V4  into  K**  V'J  +  V^,  with 


a/ 


-V--. 
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Vfd  .a  )  ■  J  )  ♦  (34«) 

and 

where  we  set  f^R,Soi  and  i‘~{6k/2k)R,S9t.  For  the  monochromatic  case  ( 
Sk  =0  ).  Rumsey  showed  that,  over  the  important  region  of  integration,  the  f 
dependence  of  Vf  is  negligible  compared  to  that  of  V^,  giving  the  intensity  spec¬ 
trum  as  the  Fourier  transform  of  ezp[~hilEs/7^(^')].  This  is  because  the  quantity 
sq/c  is  of  the  order  of  unity  sdiile  the  important  region  erf  integration  is  smaller 
than  the  scattering  disk.  When  6k  is  not  zero  (  but  small),  the  high  frequency 
behaviour  of  the  cross-spectrum  is  still  controlled  by  Vf  for  the  same  reasons. 
It  is  important  to  realize  that  these  terms  do  not  dominate  the  behaviour  of  V4 
for  all  and  £,  but  «ily  for  Unlike  the  monochromatic  case, 

when  ft  is  zero  V4  is  non-zero.  However  V4  still  increases  with  increasing  fi",  so 
that  the  Fourier  transform  integrand,  is  largest  when  /F  is  small  The 

approximate  domain  of  in  which  the  Fourier  integrand  is  large,  is  where 
In  this  regime  if  is  roughly  independent  ot  fi",  so  that  the  integrand  is 
dominated  by  the  V^.  The  limiting  value,  p'^so,  leads  to  the  requirement 

(35) 

wherein  we  expect  the  fl'  behaviour  of  t  to  be  controlled  by  the 

terms.  This  is  verified  explicitly  in  Appendix  1. 

We  proceed  in  our  analysis  by  rewriting  V4  =  1^  -f  where 


Vf®(f  J')  =  >f0-=o) 


and 


(36a) 


(38b) 


The  exact  expresaion  for  the  crou*spectrum  ii  then 


^i{i.R;kx,kj)sfxp 


y-,i«r 


(37) 


exp 


But  Vf*  is  small  compared  to  V}  allowing  us  to  expand  the  final  exponential  in  a 
Taylor  series  jdelding  a  corresponding  series  for  the  intensity  cross-spectrum 


*i(jt.R',k.6k)-  £  *f^(t,R.k.6k) 


(38) 


Truncating  this  series  after  a  few  terms  we  wiU  call  the  high  frequency  approxi¬ 
mation  since  in  the  monochromatic  case  the  corresponding  series  converges 
quickly  for  strong  scattering  and  high  spatial  frequencies.  Tbe  leading  term  is 


(39) 


exp 


[^m(3  ’  +  »o*)  ♦  /?m(?  ’  -  g-/?i  *0*)] 


6k 


d*jr 


where  Vf  ®  is  explicitly 


ff’=kf  DJ^R,  eo«(l  ^  H^)  +*l  soi{l-^)-2k^iD^RnS9t)  (40) 


We  are  now  in  a  position  to  compare  this  leading  term  with  the  "fully 
saturated”  approximation,  (31).  Since  the  phase  structure  functions  saturates 
to  for  scales  larger  than  the  outer  scale  (38)  gives  the  same  result  as  (31) 
only  for  spatial  frequencies  such  that  K>{1^^/  /?>)«ff  * .  The  only  non-n^ligible 
portion  of  the  spectrum  is  for  c  <  I/Sq,  leading  to  the  conclusion  that  the 
saturated  approximation  is  only  valid  for  a  significant  portion  of  the  spectrum 
when  the  outer  scale,  is  small  compared  to  the  size  of  the  scattering  disk. 


However,  in  virtually  every  case  of  practical  intereat,  the  acatterinf  diak  ia  amolf 
relative  to  any  eatimatea  of  an  outer  acale.  In  theae  eaaea,  the  aaturated 
approximation  becomea  va'ad  only  after  the  croaa-apectrum  has  dropped  to  a 
negligible  value. 

A  common  approach  for  dealing  with  the  intenaity  croaa-apectrum  haa 
been  to  uae  (31)  but  drop  the  phaae  ahift  factor,  t'^*^^)***. 

Ihis  ezpreaaion  can  be  compared  with  (39)  and  ia  roughly  equivalent  for  a 
aigniflcant  portion  of  the  apectrum  when  is  approximately  zero.  Uaing  the 
largest  significant  k  value  (  Sq)  leads  to  the  requirement 

*f  D^Ji,  (H-^)  +  fc|  D^iR,  (l-^)-ab,*s/>^ii;)«l  (42) 

For  a  power  law  structure  function  this  expression  can  be  simplified  by  perform¬ 
ing  a  small  ik  expansion  (  valid  for  6k  <Slk  )  leading  to 

^(^)«1  (43) 

These  points  may  be  summarized  by  sa3nng  that  the  fully  saturated 
aq>proximation  should  never  be  used  for  estimates  of  the  intensity  decorrelation 
and  (41)  should  only  be  used  when  {6k/2k)D(Ii,)<Cl.  In  fact,  since  using  (39)  is 
as  easy  as  using  (41).  the  leading  order  approximation  (39)  is  the  high  spatial 
frequency  approximation  of  chmce. 

We  now  consider  the  next-to-leading  term  in  (38).  Writing  the  phase 
structure  function  (for  the  mean  wave  number)  in  terms  of  the  phase  spectrum 

/?(•)  *  27#*(*')[i  (44) 

with  fs  as  the  power  spectrum  of  phase  shift  through  the  screea  the  secrnid 


term  in  th«  series  expansion  of  exp[-^l^* ]  is 

-|-vf  •  *  (4S) 

This  yields  the  first  correction  to  the  leading  approximation,  (39), 

=  -4^e-*^/7*,(*)e‘*‘^'*^-**)  (48) 

exp  +  sin*[^^d*p-d*«* 

where  is  given  in  (40). 

We  now  turn  our  attention  to  the  low  spatial  frequencies.  When  «dCl/  i?,, 
the  jl*  Fourier  transform  has  substantial  contributions  from  everjrwhere  in  the 
integration  plane,  and  over  most  of  this  region.  P'»R,m^k.  When  this  is  the 
case,  it  can  be  seen  from  (32),  that  V’4  will  be  dominated  by  the  first  two  terms 
on  the  right  hand  side  of  that  equation  since  the  R  dependent  terms  nearly  can* 
ceL  UnUke  the  high  frequency  approximation,  however,  the  dominant  behaviour 
is  independent,  which  is  the  source  of  the  delta  function  in  the  saturated 
approximatitm,  (31).  By  constriiction.  </,>=</9>=l  since  we  set  the  input 
amplitude  to  unity.  Since  the  dependent  terms  in  are  small,  the  eiq>onent 
may  be  eiq>anded  in  a  Taylor  series 

*)=exp  (47) 

yielding  the  low  frequency  series  for  the  cross-spectrum 

♦/(«,/?;* ,.*a)=  t 

nmo 


(48) 


Insertinf  the  above  expansion  into  (31)  yields  explicit  expressions  for  the 
temu  in  the  low-frequency  series.  Using  (44)  to  express  the  structure  functions 
of  (47)  in  terms  of  the  phase  spectrum,  we  find  that  the  two  leading  terms  in  the 
low-frequency  series  ve 

♦/!?>(*./?:iki.*s)*«(*)  (♦») 


and 


4t(g)|8in*(/^  So**/  2)-sin*(^/^  *o**/4)j 


(50) 


I 


As  mentioned  earlier,  the  delta  functitm  is  the  same  as  in  the  saturated  approxi¬ 
mation.  but  the  behaviour  for  0<*«l/i^  is  very  different 
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In  this  appendix  we  show  that  ^  dependence  of  Vf  is  negligible  compared 
to  that  of  Vj  where,  from  (27). 

(Ala) 

and 

(|’.J'.t>=*  fD^if*^‘)*kiD^f-i‘)^k  ik^Djt^*f)*D^-f)]  (Alb) 

We  have  argued  in  section  4  that  the  dominant  region  of  integration  is  for  small 
^OCi?a  sqc.  In  that  case,  the  two  terms  in  the  brackets  in  (Alb)  may  be  expanded 
as 

So*-^)=2/?^(^,  So«)+ 

The  second  term  may  be  bounded  by 

l>-if  DJii)  la.i«;,,.]=^^l  *o«)  I  (A3) 

giving  the  leading  ^  dependence  of  1^  as 

if-*  *0*)  (A4) 

In  verifying  that  Vj  controls  the  ^  dependence  we  will  examine  three  cases 

Qise  1:  Sqk<CB 

Qise2:/r<C^^So«  (A5) 

Chse  3:  *o**^ 

In  the  first  case,  fi  is  large  compared  to  M/i^Soc/2fc  giving  the  leading 
dependence  of  as 


M 


(A2) 


K?~2*,*3/?m(^) 


(A6) 


103 


to  that  we  require 

(A7) 

Using  a  model  power-law  structure  function  as  in  (30)  this  requirement  becomes 


iJf’M 


-z  2 

«  * 


p(p-l) 


(A8) 


Snce  p  is  t3rpicaUy  in  the  range  Q<p<2,  we  find  an  upper  bound  for  the  left-hand 
side  by  setting  fi'-Sc  and  c=Sq  ‘  resulting  in  the  requirement 


(A9) 


Therefore,  the  term  controls  the  dependence  when  the  strength  of  scatter¬ 
ing  parameter,  D(R  sub  s  ).  is  large.  When  p  is  1  or  2.  our  argument  breaks 
down.  For  p=l  the  approximation  still  holds  by  another  argument,  while  for 
p=2.  there  are  no  scintillattions. 

In  the  second  case,  fi'  is  small  compared  to  6kR,s^K/^,  so  that  a  Taylor 
expansion  of  the  term  yields  the  leading  /T  dependence 

,6k 


Zk 


(AlO) 


This  leads  to  the  requirement 


*o«) 


which  for  a  power-law  structure  function  is  equivalent  to 


(All) 


1<C 


2fc 


(A12) 


Since  p  is  less  than  2.  this  implies  that  6k  is  small  compared  to  tvnce  the  mean 
of  the  wavenumbers. 


Finally,  in  the  third  case,  f  is  about  the  same  size  as  6kR,  Sqc/  2k  so  that 


the  requirement  for  neflecUng  the  ^  dependence  of  Vf  becomei 

*  Ak 

(A13) 

which,  for  a  power*law  structure  function  leads  to 

(AM) 

and  is  easily  satisAed  for  small  wavenumber  differences. 


In  all  three  case  considered  above,  the  combination  of  structure  func¬ 
tions  in  V9  dominate  the  dependence  if  D(ffaau^  s  )  is  large  and  6k/2£  a 


The  Scintillation  Theory  of  Eclipse  Shadow  Bands 

by 

Johanan  L  Codraia  * 


ikatnet.  The  resuha  of  a  tbeoretieal  iimatifatioa  of  aolar  cchpae  ahadow  handa  are  preaenUd.  The 
atiidy  providea  both  quantitative  and  qualitative  inaicht  into  the  factor*  f overning  the  viaibility  of 
ahadow  banda.  Without  uae  of  auxiliary  l^potheaea  aueh  aa  aolar  hmb  darkeniiv  or  tur¬ 

bulence  enhancementa  cauaed  by  the  paaaage  of  the  hmar  ahadcnr,  aD  of  the  aalient  featurea  of  the 
ahadow  banda  arc  explained  by  atandard  acintillatiaD  theoey.  The  oontraal  ia  found  to  be  greater  for 
ahorter  wavelengtha  and  the  band  apacing  to  acale  like  the  aquareroot  of  the  wavelength  very  near 
totality.  For  timea  greater  than  about  10  aecooda  before  (or  after)  totality  the  band  apacing  becomea 
frequency  independent  and  the  scintillations  are  dominated  by  turbulence  near  the  ground.  The  tur- 
btdenoe  atai^  reip<Hiaib}e  fmr  ahadow  banda  ia  found  to  be  below  two  Idloineters  in  altitude.  Tur¬ 
bulence  at  the  tropopause  is  found  to  have  no  intact  on  shadow  bands  until  2-3  seconds  from  total¬ 
ity.  Longer  eclipses  are  expected  to  show  bands  with  greater  contrast  and  linearity.  Intensity  corre¬ 
lation  scales  are  typically  leas  than  10  cm  within  30  seconds  of  totality.  The  scintillation  theory  pred¬ 
ictions  for  ahadow  band  structure,  motion,  and  evolution  are  found  to  be  in  renwrkable  agreement 
with  both  visual  and  photoelectric  observations. 

1.  Introdnetioa. 

In  the  minutes  just  preceding  and  just  following  totality  during  solar 
eclipses,  faint  shadowy  patterns  have  sometimes  been  seen  to  appear  projected 
on  the  ground  [l].  Hiese  "shadow  bands"  are  initially  seen  to  be  random  and 
wavy,  becoming  more  linear  and  aligned  with  the  edge  of  the  lunar  shadow  as 
totality  approaches.  At  the  same  time  the  band  spacing  decreases  and  the 
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contrast  increases.  The  bands  appear  to  move  in  a  drection  perpendicular  to 
their  elo.\2ation  [l].  Photoelectric  studies  have  shown  that  the  bands'  contrast 
is  greater  at  shorter  wavelengths  while,  near  totality,  the  typical  width  of  the 
bands  increases  with  wavelength  [2].  The  phenomeni  appear  to  reverse  their 
general  trends  after  totality  with  the  exception  of  the  band  motion,  which 
appears  to  be  unrelated  to  the  locatim  of  the  lunar  shalow. 

It  has  long  been  suggested  that  shadow  bands  are  most  likely  a  scintillation 
phenomenon  caused  by  the  Earth's  atmosphere  [3].  However,  details  such  as 
whether  or  not  shadow  bands  will  be  seen  and  how  wdl,  the  wavelength  depen¬ 
dence  of  the  bands,  and  the  nature  ol  the  atmospbeac  irregulariUes  have  not 
been  well  understood.  The  wavelength  dependence  hat  been  attributed  to  stdar 
limb  darkening  [4]  and  the  turbulence  to  possiUy  have  been  enhanced  by  the 
thermal  effects  of  the  moon's  shadow  [5].  In  this  paper  it  will  be  shown  that  all 
of  the  effects  of  shadow  bands  can  be  explained  with  standard  scintillation 
theory,  without  the  use  of  limb  darkening  (m*  enhanced  turbulence. 

In  the  next  section,  the  problem  of  a  general,  extended,  incoherent  light 
source  at  infinity  illuminating  a  thin  phase  screen  wil  be  considered.  A  model 
for  the  two  point  intensity  correlation  will  be  found  is  well  as  the  normalized 
intensity  variance  (the  scintillation  index).  In  section  3  the  speci&c  case  of  cres¬ 
cent  source  distributions  will  be  presented.  Section  4  sill  utilize  the  thin  screen 
model  to  find  the  shadow  bands'  intensity  variance  ant  wavelength  dependence. 
In  section  5  the  morphology  and  scales  of  the  thin  screen  shadow  bonds  will  be 
discussed.  Ihe  intensity  fluctuations*  dependence  on  time  and  wavelength  will 
be  calculated.  In  section  6  the  theory  will  be  generalked  to  consider  extended 
atmospheric  turbulence.  The  implications  of  the  eAended  atmospheric  tur¬ 
bulence  model  for  the  shadow  bands’  wavelength  dependence,  contrast,  and 
structure  will  be  showa  Section  7  will  summarize  the  results  of  the  theory  and 
make  comparisons  to  data. 

2.  Ihe  niaae  Screen  Model 

The  shadows  on  the  ground  can  be  described  as  a  random  intensity  pattern, 
1(1),  at  any  given  moment.  Here,  1  is  a  position  vector  in  the  observation  plane. 


The  inteniity  patterns  result  from  ten^rature  and  density  fluctuations  in  the 
atmosphere  and  change  in  a  semi-random  fashion  as  the  irregularities  evolve 
and  are  convected  hy  the  wind.  We  will  treat  the  light  as  being  essentially  mono¬ 
chromatic.  This  will  be  a  good  approximation  if  the  intensity  is  detected 
through  narrow-band  Alters.  A  study  of  the  two-point  intensity  correlation, 
Q(fl)  -  ***  fl)>.  will  lend  a  great  deal  of  insight  to  the  nature  of  the  sha¬ 

dow  bands.  The  angle  brackets  denote  ensemble  averaging  over  fluctuations  in 
the  atmosphere.  The  intensity  correlation  function  allows  us  to  discuss  how  the 
correlation  scales  vary  in  time  as  well  as  the  variance  in  the  intensity.  The  vari¬ 
ance  is  <(/-</>)*>  =  C/(0)  -  C/(«b).  The  "scintillation  index”  is  the  normalized 
rms  fluctuations  in  intensity,  m  =  <(/-</>)*>*/  <I>.  Therefore, 


_  C/(o)-g/H 

■  Q(-) 


(2.1) 


In  this  section  we  derive  C/(t)  for  a  weakly  scattering  layer  of  turbulent  atmo¬ 
sphere  at  height  z  (a  thin  "phase  changing  screen”).  Below  this  layer  the  index 
of  refraction  is  constant  Without  loss  of  generality,  the  sun  vnU  be  considered 
to  be  at  the  zenith  (flgure[]).  The  index  of  refraction  fluctuations,  n(r),  are 
assumed  to  arise  from  a  turbulent  process  in  the  air.  The  "turbulence  spec¬ 
trum”  is  defined  as 


♦n(fl)  =  (2.2) 

where  is  the  three  dimensional  wave  vector  and  the  index  of  refraction  autoco¬ 
variance,  B^(f)  =  <n(flo)n(fo  >  f)>,  is  assumed  to  be  hcmiogeneous  and  isotropic 
within  the  turbulent  region.  For  spatial  frequencies  important  to  scintillations, 
we  can  treat  the  turbulence  spectrum  as  a  power-law  with  spectral  index  a  [6]. 

♦•.(4)  =  0.033qfg-«  (2.3) 

where  €$  is  the  turbulence  structure  constant  and  a  =  11/3  fn*  the  case  of  Kol¬ 
mogorov  turbulence  [6]. 

Propagation  of  light  through  the  turbulent  layer  and  to  the  ground  can  be 
treated  as  narrow  angle  forward  scattering.  In  situations  where  the  radiation  is 
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flcaiterad  into  •  norrow  cone  About  tho  incoming  direction  and  the  fluctuations 
in  the  index  of  refraction  are  much  leaa  than  unity,  the  Helmholtz  wave  equation 
reduces  to  the  "parabolic  wave  equation" 

ff-=  ^9‘E-ikni%M)E  (ZA) 

where  B-Re(£’e*t^ is  the  scalar  electric  fleld.  z  is  the  coordinate  in  the 
direction  of  propagation,  i  is  a  eecUx*  perpendicular  to  the  z  axis,  k  sZn/\  is 
the  wavenumber,  X  is  the  wavelength,  and  is  the 

transverse  Laplacian  [6].  In  the  region  between  the  screen  and  the  ground 
n  s  0.  Setting  n  s  0  in  (2.4)  and  Fourier  transforming  with  respect  to  t  yields 

«5) 

Os  aT®  '  ' 

where 


(2-6a) 


Ei%z)  =  7f(«.s)s«**d*« 


(2.6b) 


and  t  is  the  two  dimensional  spatial  frequency.  Equation  (2.5)  implies 


where  S{tt,0)  is  the  Fourier  transform  of  the  fleld  just  under  the  turbulent  layer 
(s  s  0^).  This  allows  us  to  write  the  fleld  on  the  ground  as 


E{Xz)  = 


*.l-  d*« 

2k 


Define  Eo(t)  to  be  the  fleld  at  s  s  0*^  wMch  would  result  from  a  plane  wave  imp* 
inging  on  the  phase  screen  along  the  s  axis.  Since  the  screen  is  assumed  to  be 
thin,  a  wave  that  is  incident  on  the  screen  with  an  angle  4  to  the  s  axis  will 
result  in  a  fleld  at  s  =0*  of  £*0(1)  exp(Hib^*t).  Using  equation  (2.7)  and  the 
Fourier  shift  theorem  we  write  the  fleld  on  the  ground  resulting  from  a  unit 


"T  *_  "*•  *^1*  **■  *’*"*_*^  *1  •  ".."i 
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amplitude  plane  wave  incident  from  an  an(le  4  aa 

E{±m:4)  =  fioii‘k4)  e*p|l*.t .  (l««  (^8) 

The  corresponding  intensity  pattern  is 

sdiere  the  asterisk  (*)  denotes  complex  coQjugatioiL  Using  equation  (2.S).  this  is 
written  in  integral  form  as 

-  ff  exp|it.(tt  -  <t)  -  -  *f)j  (2*9) 

-  k4)g^kt  -  k4)  d*ci  d*«a 

If  we  illuminate  the  screen  by  an  incoherent,  extended  source  with  a  bright* 
ness  distribution.  B(4),  the  intensity  pattern  becomes 

/(t*)  =  f  B{4)I(%*:4)  (2.10) 

^ce  Bi4)  is  zero  beyond  the  edge  of  the  source,  the  limits  of  integration  can 
be  written  as  infinite  without  any  difficulty.  An  examination  of  equation  (2.9) 
shows  that 

I{Xx;4)  -  /(t+s^,s.-^  =  0)  »/^(l+ s^)  (2.11) 

where  the  subscript,  pt,  refers  to  a  point  brightness  distribution  at  ^  =  0.  Equa¬ 
tion  (2. 11)  is  valid  as  long  as  the  angle  4  ia  small  compared  to  a  radian.  Combin¬ 
ing  equations  (2.10)  and  (2.11)  we  write 

/(*)  =  7^(^)/^(l+s^)d*d  (2.12) 

leading  to  the  intensity  correlation 

C/(i)  =  ff  B{4x)B{4^c:^{t-t{4y  -4,))  d«4,  d^,  (2.i3) 

where  f^(l)  -  </yt(f)/^(t -^  •)>  and  translational  invariance  has  been 


•itumed.  Note  that  C/(f)  need  not  be  ieotropic  if  B{^)  ii  anisotropic.  This  is 
why  the  scintillations  of  stars  and  planets  are  isotropic  [7]  but  the  shadow  bands 
are  not  The  intensity  spectrum  is  defined  by 

♦/(*)  =  ^2^/C/(i)e-«-*d«s  (2.14) 

Therefore,  equation  (2. 13)  leads  to 

♦,(«)  =  (2.15) 

where  ^(sd)  is  the  Fourier  transform  (rf  the  brightness  distribution  with  respect 
to^  and  *^(t)  is  the  Fourier  transform  of  C^(l).  Tlus  is  the  Cohen-Salpeter  for¬ 
mula  for  the  intensity  spectrum  of  an  extended,  incoherent  source  [B]. 

The  intensity  spectrum  in  equation  (2.15)  consiats  of  two  factors.  There  is 
the  intensity  spectrum  of  a  point  source.  and  the  source  hinction. 

tff(s£)|* .  which  smoothes  the  intensity  fiuctuations.  of  an  inten¬ 

sity  spectrum  wUl  reveal  a  delta  funcUon  at  the  origin  due  to  the  mean  intensity. 
The  remainder  of  the  q>ectrum  is  the  Fourier  transform  of  the  intensity  autoco- 
variance  and  is  a  consequence  of  the  intensity  fiuctuations.  Writii^  the  fiuctua¬ 
tions  in  intenaty  as  /'  =  f  -  </>  we  note  that 

♦/(£)  =  </>*d(*)  +  ♦/.(£)  (2.16) 

where  d(£)  is  the  two  dimensional  Dirac  delta  function  and  #/■(£)  is  the  intensity 
fluctuatim  spectrum.  Combining  equations  (2.1)  and  (2.18).  the  scintillation 
index  becomes 


FVom  equaUons  (2. 15)  and  (2. 17)  we  can  see  why  the  scintillation  strength  varies 
with  source  size.  A  brightness  distribution  with  a  characteristic  width  Bq  will 
have  a  Fourier  transform  with  a  characteristic  width  of  order  2ir/do .  TTierefore. 
the  brightness  distribution  acts  as  a  low-pass  filter  in  c  with  a  cutoff  around 
2n/(sdo)  .  As  we  increase  the  size  of  the  source,  we  decrease  the  cutoff 


c 


^quvncy  in  4/.  Thii  diminishei  the  integral  in  equation  (2.17)  leading  to  a 
decreaie  in  the  icintillation  index,  m*.  ITut  is  the  reason  that  stare  twinkle  so 
much  more  than  planets.  The  situation  for  shadow  bands  is  complicated  by  the 
fact  that  the  brightness  distribution  for  a  solar  crescent  is  so  as)rmmetrical.  In 
such  a  case,  the  Altering  becomes  a  fimction  of  direction.  This  will  be  discussed 
in  detail  in  the  following  sections. 

In  addition  to  the  effect  of  the  source  on  the  scintillations,  we  must  also 
consider  the  point  source  intensity  spectrum.  This  function  describes 

how  a  star  would  scintillate  under  the  same  conditions.  To  discuss  this  function, 
it  is  necessary  to  look  at  the  fourth  moment  of  the  electric  Aeld.  Tlus  approach 
is  necessary  because  there  is  no  completely  correct  propagation  equation  for 
the  intensity  correlation  function.  Die  fourth  moment  is  deAned 

where  the  E{%^)  factors  are  the  electric  Aeld  at  the  points  resulting  from  a 
plane  wave  impinging  perpendicudarly  on  the  screen.  The  intensity  correlation. 
qe(i).  is  found  from  the  fourth  moment  as 


s  17(1.1.1 +1.1  +  1) 


(2.1B) 


The  parabdic  wave  equation,  equation  (2.4),  can  be  used  to  derive  a  propagation 
equation  tor  the  fourth  moment.  If  [6].  In  the  region  below  the  screen,  where 
n  =  0.  this  equation  is 


^  -  V|  +  V|  -  7D^(ti.V«s.t^*) 


(2.19) 


where  VjJ  is  the  transverse  Laplacian  with  respect  to  Following  Rumsey  [9], 
we  change  variables  to  A,  f,  and  B  acc«x‘ding  to  the  transformation 
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1-1  1-1 


(2.20) 


^ce  both  the  medium's  statistical  properties  and  the  geometry  of  the  problem 
are  Invariant  under  translations  in  1,  we  can  assume  that  W  will  be  independent 


of  A.  Tbut.  oquotion  (2. 19)  bocoinot 


(z-zt) 


where  we  have  set  A  s  J  s  0  without  any  toss  of  generality.  This  leaves  the  four 
points  in  a  parallelogram  with  sides  $  and  %  The  intensity  correlation  will  be 
found  by  setting  f-Q.  For  convenience,  a  superscript  $  or  f  will  denote  a 
Fourier  transform  with  respect  to  that  variable.  The  CMresponding  spatial  bw- 
quency  variable  will  be  written  as  a  subscripted  A.  Transforming  equation  (2.21) 


we  get 


g>  ^ 


This  implies  that 


ir»(/^z,«  =o‘) 


Inverse  Fourier  transforming  with  respect  to  we  find 


9^(A.f.s)  * 


(2.22) 


where  we  have  dropped  the  P  subscript  on  A  and  the  subscript  0  on  refers  to 
s=0*. 

To  find  Wq,  we  assume  that  the  screen  contributes  a  random  phase  shift. 
0 

at  the  point  %  where  I  is  the  thickness  of  the  screen. 

Snce  we  are  illuminating  the  screen  with  a  normally  incident  plane  wave,  the 
electric  field  at  t  =0*^  would  be  unity  in  the  absence  of  the  screen.  With  the 
screen,  the  field  is  exp[iF(A)].  Therefore,  a  realiation  for  Ag  is  ezp[i(F(^)  • 
F(^  f(^)  '  ^(At))].  Assuming  f(ti  to  be  a  aro  mean  Gaussian  random  vari¬ 

able.  we  write 


lfo(^,Ai.A}.ti)  *  eap  -  ^f^is  +  ^  *  ^is  *  (2.23) 


where  ~  <(f(^)  *f(%))*>  u  the  phaw  structure  function  of  the  screen  [6]. 


Chanfini  variabiM  in  aquation  (&23)  according  to  aquation  (2.20)  with  d s I  >0. 

wa  gat 

H'oili)  =  exp[-  V(0)  -  D(f)  >  *f)*  D(^  ■  f))]  (2.24) 

Examination  of  the  exponent  in  equation  (2.24)  shows  that  it  is  bounded  between 
-D(f)  and  0.  In  weak  scattering,  the  exp<mential  of  the  structure  functions  that 
depend  on  ^  in  equation  (2.24)  can  be  Taylor  expanded  to  give 

-  f^»)[l.X?0)+Jf(i>(?  +  f)4/)0.f))4  ) 

The  exp[  -  X>(f)]  factor  may  be  ignored  in  weak  scattering.  Assimung  that  the 
screen  is  thicker  than  several  index  of  refraction  correlation  scales,  we  write  the 
phase  structure  functicxi  as 

/7(1)  s  4rnkHf^(t.q»  =0)[l  •  t«**ld*«  (2.25) 

where  q,  is  the  vertical  component  of  the  wave  vector  4  Substituting  equation 
(2.25)  in  the  expansion  above,  we  find  after  some  manipulation  that 

"'^(d.f)  =  «(*>  8»rl;*fsin»[^*,(«) 

Combining  this  result  with  equations  (2.22)  and  (2.18)  we  see  that 

♦^(2)  =  «(«)  +  8Tr**lsiii*[^*,(*)  (2.26) 

Hus  is  the  weak  scattering  iq}proximation  to  the  intensity  spectrum  for  a  point 
source.  The  sin^(sx^/2fe)  facUxr  is  called  the  "fVesnel  filter."  The  general  inten* 
sity  spectrum  is  found  by  cmnbining  equations  (2.26)  and  (2.15)  to  give 

♦/(*)  =  I&uu6{£)  +  {Zn)*bnkH  l5'(s2)|*  sin*[^*»(*)  (2.27) 

where  we  have  made  use  of  the  fact  that  (2n)*^(0)  s  /mm.  the  total  intensity  of 
the  extended  source.  The  intensity  fluctuation  spectrum  is  the  second  term  in 
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W«  can  now  write  down  the  intenaity  eorroletion  on  the  ground  and  the  scin* 
tiUation  index  IVom  equationi  (2.27)  and  (2. 18)  we  make  the  identification 

</>  = 


Therefore  the  intensity  correlation  is 


Ciit)  =  /&-  +  (2ff)*8nls«l7e  -<«**|fi'(s8)|*  sin* 


^♦,(*)d*c  (2.28) 


and  the  scintillaUon  index  is 


■sm' 


(2.29) 


EquaUons  (2.28)  and  (2.29)  are  valid  for  essentially  any  broad  source  seen 
through  a  tt^  turbulent  layer.  The  crescent  source  distribution  is  what  makes 
the  special  case  of  shadow  bands  so  umque.  We  look  at  this  next 


3.  Orescent  Source  Diatrihutiaiks. 

The  very  thin  solar  crescent  is  i^at  makes  the  common  scintillation  effects 
in  shadow  bands  so  remarkable.  Its  sheer  brilliance  sets  it  apart  from  other 
scintillating  objects  in  that  it  may  be  (4>served  by  backscattering  off  the  ground. 
The  crescent  also  changes  its  shape  and  size  rapidly  as  totality  nears.  The 
manner  and  rate  at  which  this  evolution  takes  place  is  governed  by  the  relative 
sizes  of  the  sun  and  moon  which  change  frmn  one  eclipse  to  the  next.  The  cres¬ 
cent  may  even  change  its  complete  nature,  as  in  annular  ecBpses,  from  a  cres¬ 
cent  to  a  ring.  During  some  total  ecl4>ses,  the  moon  may  be  very  neaiiy  the 
same  apparent  diameter  as  the  sun.  On  other  occasions,  the  moon  may  appear 
considerably  larger  than  the  sun.  So  that  we  may  study  the  relatimships 
between  shadow  bands  with  different  eclipse  geometries,  we  need  to  define  some 
characterizing  parameters.  Let  the  angular  radii  of  the  sun  and  morni  be  Rs 
and  Rg,  respectively.  Even  though  both  of  these  radii  may  change  from  eclipse 
to  eclipse,  it  is  their  ratio  which  defines  the  nature  of  the  eclipse  and  has  the 
greatest  effect  on  shadow  bands.  Therefore,  we  define  an  'eclipse  parameter’,  c. 
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where  the  subscripted  t  is  the  time  at  the  appropriate  contact  llie  parameter 
e  is  positive  for  total  eclipses  and  negative  for  annular  eclipses.  A  total  eclipse 
where  Baily's  beads  are  prominent  and  seen  over  a  large  arc  of  the  crescent  will 
have  an  t  very  near  zera  A  long  total  eclipse  will  have  a  larger  value  of  e. 
Another  important  parameter  is  the  time  away  from  totality.  Define  the  angular 
velocity  of  the  moon  relative  to  the  sun  as  (I  The  normalized  time  away  from 
totality  is  defined 

^  ^3^2) 

/fs 

where  fcmiaei  is  time  of  either  second  or  third  contact,  depending  on 
whether  we  are  in  the  waxing  or  waning  phases  of  the  eclipse.  At  both  first  and 
fourth  contact,  r  has  a  value  of  two.  At  second  or  third  contact,  r  is  zero.  As  a 
benchmark  number,  we  will  use  60  minutes  for  the  interval  f//  -  f/.  This  gives 
approximate  formulae  for  r  and  e  as 

duraticm  of  totality  in  minutes 
80 

time  away  from  totality  in  minutes 


This  gives  typical  values  for  e  as  0  to  0.12  for  total  eclipses.  Shadow  bands  are 
usually  seen  for  times  less  than  two  minutes  away  hom  totality.  Ilierefore,  the 
values  of  r  that  we  are  interested  in  lie  in  the  range  0  to  about  0.1 . 

A  set  of  solar  crescents  is  shown  in  figure  1.  The  most  striking  dillerence 
between  eclipses  with  different  values  of  e  are  extent  of  the  crescent’s  horns. 
Ti*eating  the  moon  and  sun  as  circles  whose  centers  are  separated  by  an  angle 
/?5(e  -f  r)  we  write  the  horn  angles  from  the  center  of  the  moon  and  the  sun  as 
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ThM«  angles  are  shown  in  figure  2.  Figure  3  shows  the  evolution  of  the  solar 
horn  angle,  fs,  with  r  for  verioua  valuea  of  epeilon.  Eclipaes  with  smaller  values 
of  e  tend  to  keep  their  horns  up  to  the  last  moments  before  totality  while 
eclipses  with  larger  values  of  epsilon  lose  their  horns  rapidly  and  display  a 
lengthy  "Diamond  Ring"  effect.  Smaller  crescent  born  angles  mean  that  the 
solar  crescent  is  forming  a  better  "slit."  Thus,  we  might  expect  that  eclipses 
with  larger  values  of  t  would  make  crisper  shadow  bands.  This  is  true,  but  the 
horn  angle  will  affect  the  observed  structure  of  the  bands  more  profoundly  than 
the  contrast.  Ihe  other  point  th&l  we  need  to  consider  is  the  mean  intensity  of 
the  crescent  Relative  to  the  sun’s  full  disk  intennty  we  write  the  eresent’s 
intensity  as 

0  ~  •  sin(2pir))  (3.4) 

Ihis  normalixed  intensity  is  shown  in  figure  4.  The  larger  horns  in  small  e 
eclipaes  keep  the  intensity  up  much  higher  than  in  eclipses  with  a  larger  value 
of  t. 

The  center  of  the  crescent  lies  along  the  line  between  the  apparent  centers 
of  the  sun  and  moon  Therefore,  if  an  observer  is  not  on  the  center  of  the  path 
of  totality,  the  crescent's  orientation  will  be  a  function  of  time.  Ihis  effect  is 
illustrated  in  flgwe  5.  If  the  observer  is  not  on  the  center-line,  while  s  is 
unchanged,  r  will  become  a  more  complicated  function  of  time  than  given  in 
(3.2). 

(tJ  +  2eToVi  -p*  +  e*)* - e  o^p^l 

^  (Tj  +  p*e*)»-e  P>1 

where  tq  is  the  value  of  r  that  would  be  found  by  using  (3.2)  and  p  is  a  normal¬ 
ized  "impact  parameter"  showing  the  observer's  fractional  distance  from  the 
center  of  the  path  of  totality  toward  the  edge.  If  p>l  the  observer  will  not 
experience  totality  at  all.  In  that  case.  Tq  is  calculated  using  (3.2)  with 
replaced  with  the  time  of  greatest  partial  phase.  Note  that  for  p>l.  r  never  goes 
to  zero  (figure  6). 


Calculation  of  the  inteneity  ipeetrum  requirei  knowledge  of  the  ereecent's 
angular  power  spectrum.  .  This  is  very  difficult  to  calculate  exactly. 

Fortunately,  we  can  make  two  good  approximations  to  |^(sg))‘  :  one  that  is  valid 
for  large  values  of  eg  and  one  that  is  good  for  small  values  of  sg  .  Ignoring  the 
effects  of  solar  limb  darkening,  we  can  write  the  brightness  distribution  as 

^  I  ffo  if  ^  is  inside  the  crescent 
=  I  0  otherwise 


The  corresponding  source  spectrum  is 


where  C  denotes  that  the  region  of  integration  is  over  the  crescent.  Using  a  two 
dimensional  form  of  the  divergence  theorem  we  arrive  at 


(4.1) 


where  dC  denotes  the  boundary  of  the  crescent.  n(d)  is  the  unit  outwardrfacing 
normal,  and  old  is  the  differential  angular  path  length  along  the  edge  (rf  the  cres¬ 
cent.  Equation  (4.1)  lends  itself  to  asymptotic  analysis  as  xkRs  becomes  large. 
Using  the  method  of  stationary  phase  and  the  fact  that  we  can  write 


{Zn)\xKRsf 


zkRs 

2 


<e  -  (e  +  T)cos(f )) 


(4.2) 


forsic/?5»l  and 


where  f  is  the  angle  between  g  and  the  axis  perpendicular  to  the  crescent 
(figure  2).  For  i>>^s  •  l^(*^)l^  Is  essentially  zero.  The  region  F5<‘lf'<Fjr  a 
different  asymptotic  behavior  but  may  be  neglected  in  practice.  Ifhen  we  use 
(4.2)  in  (2.29)  to  get  the  scintillation  index,  we  will  be  integrating  over  if.  If 

fs  <  tbs  result  of  this  integration  is  given  asymptotically  as 


f\B{*(trdif 

9 


ISnBiRiil  4  <)*  _  coajMKRsT  +  ^ 

(2iry«(a*^,)»  (r«*i?,(s  +  T))* 


(4.3) 


120 


for  mkFs^I-  An  approximate  form  for  for  the  crescent'a  power  spectrum  alonf 
the  if*  s  0  axis  is  given  by 


^  exp  tt3(exp 

U  S  XkRsT 

Oi  S  0.9 


-a2U*  a,u 

®s 


•  1)  -OiU 


02  =  exp(.027ln*(e  +  t)  -  1.7Bln(e  +  t)  -  2.83) 
og  =  3exp(-11.0(e  +  t)) 


(4.4) 


Equation  (4.4)  is  the  result  of  numerical  curve  fitting  and  is  good  for  xkRs^ 
between  zero  and  about  3  or  4.  Numerical  calculations  show  that  (4.2)  and  (4.4) 
form  an  approximation  that  is  good  tor  the  entire  range  of  c  fcr  ^  =  0  (figure  7). 
The  integration  will  be  approximated  by  multiplying  (4.4)  by  the  factor  4^5. 

We  now  have  all  of  the  information  necessary  to  find  the  scintillation  index 
due  to  a  thin  screen.  We  need  only  to  pick  the  wavelength,  A.  the  distance  to  the 
screen,  z ,  the  strength  of  turbulence,  C^.  and  the  parameters  e  and  r.  To  see 
which  altitudes  are  most  efficient  at  generating  intensity  fluctuations,  we  define 
a  "scinUllation  efficiency" 

T7^t(*.A.e,T)  js  (4.5) 


where  is  the  result  of  integrating  (2.29).  larger  values  of  imply  that 
more  intensity  variance  will  be  gernerated  for  a  given  level  of  turbulence,  figure 
8  shows  vs.  z  for  different  values  of  X,  e,  and  r.  At  larger  values  of  s,  the 
source  spectrum  limits  the  integral  in  (2.29).  This  "source  averaging"  reduces 
the  scintillations  by  smearing  the  pattern  on  the  ground.  When  the  source  width 
is  much  broader  than  VA/  s  .  the  sine-squared  factor  may  be  approximated  as 
(z/c^/2k)^.  This  provides  a  factor  which  cancels  the  in  front  of  the  integraL 
Therefore,  when  the  turbulence  is  sufficiently  distant  from  the  (^server  or  the 
source  sufficiently  broad,  the  scintillations  will  not  have  a  wavelength  depen¬ 
dence.  At  lower  altitudes,  vdicn  the  source  does  not  provide  a  cutoff,  the  inten¬ 
sity  variance  will  be  greater  for  shorter  wavelengths.  This  is  because  the  Fresnel 
frequency  is  proportional  to  VX  and  the  Fresnel  filter  gives  a  larger  contribution 
to  (2.29)  for  shorter  wavelengths.  The  transition  altitude  between  wavelength 


dependent  end  wevelength  independent  scintilletions  dependi  on  the  width  of 
the  crescent  end  is  proportional  to  X/  r*.  The  eclipse  parameter,  s.  also  has  an 
effect  on  the  scintillation  efficiency.  While  not  as  important  as  source  averaging, 
a  smaller  value  of  e  will  degrade  the  scintillation  efficiency  at  all  altitudes.  For 
two  eclipses  with  identical  viewing  situations  and  weather,  an  eclipse  with  es  .01 
will  have  rms  intensity  fluctuations  that  are  20%  less  than  an  eclipse  with  e  = .  10. 
Because  shadow  bands  are  usually  near  the  limit  of  the  eye's  ability  to  detect 
contrast  and  eclipses  with  smaller  values  of  e  are  brighter,  the  e  dependence 
coidd  mean  the  difference  between  whether  the  bands  are  seen  or  not. 

The  scintillation  efficiency  of  a  star  is  never  limited  by  its  source  distribu¬ 
tion.  It  is  given  by 

which  increases  monotonically  with  z.  It  is  because  distant  turbulence  is  more 
effective  than  nearby  turbulence  that  stellar  scintillations  are  commonly  said  to 
arise  in  the  tropopause  at  an  altitude  of  10  -  12  km.  However,  lower  altitude  tur¬ 
bulence,  which  is  more  responsible  for  the  related  phenomenon  of  "seeing,"  can 
contribute  substantially  to  a  sUir's  scintillations.  Young  [7]  noted  that  plane¬ 
tary  scintillation  intensity  spectra  could  not  be  well  modeled  by  considering  only 
a  thin  screen  at  the  tropopause.  A  finite  diauneter  source  may  be  more  sensitive 
to  high  level  or  low  level  turbulence  or  both  depending  on  the  scintillation 
efficiency  profile.  Shadow  bands  are  unique  in  that  their  scintillation  efficiency 
evolves  rapidly  and  dramatically  with  time.  For  larger  values  of  r,  shadow  bands 
are  more  sensitive  to  turbulence  near  the  ground.  As  totality  nears,  turbulence 
at  higher  altitudes  begins  to  play  a  more  important  ri^e.  Hus  will  change  the 
shadow  bands'  wavelength  dependence  as  well  as  their  morphology  as  time 
progresses.  This  behavior  is  included  with  a  model  for  the  distributed  atmos¬ 
pheric  turbulence  in  section  6. 

5.  Riase  Screoi  Shadow  Band  Structure. 

Features  of  shadow  bands  such  as  their  orientation,  spacing,  motion,  and 
linearity  are  all  able  to  be  discussed  in  terms  of  the  spatial  intensity  fluctuation 


spectrum  (2.27).  This  spectrum  is  en  ensemble  everege  of  the  instenteneous 
spetiel  spectre  end  is  ergodlcelly  releted  to  e  spetial  everege.  nef  e  temporel 
everege.  Since  the  spectrum  evolves  very  rapidly  near  totality,  temporal  statis* 
tics  are  at  best  confused  and  at  worst  meaningless.  However,  the  evolving  mor* 
phology  of  the  bands  may  be  understood  as  snap-shot  realizations  of  the  spatial 
spectrum  calculated  tor  the  appropriate  time.  The  simplest  such  insight  is  the 
orientation  of  the  bands  with  time.  As  mentioned  in  section  3.  the  crescent's 
orientation  is  a  function  of  time  if  the  observer  is  not  on  the  center  of  the  path 
of  totality.  The  source  power  spectrum.  |^(sg)|*.  rotates  with  the  source. 
Iberetore.  since  |ff(sg)|*  is  the  «ily  source  of  anisotropy  in  #/(2).  the  shadow 
bands  will  rotate  with  the  source.  A  posable  caveat  is  when  there  is  some  aniso¬ 
tropy  in  the  turbulence  spectrum  not  included  in  the  standard  Kolmogorov  tur¬ 
bulence  model  Such  an  anisotropy  might  be  due  to  atmo8|dieiic  gravity  waves 
propagating,  through  the  background  tuii}ulence.  Although  the  medium's  aniso¬ 
tropy  is  a  possibility,  it  is  unlikely  to  be  an  in^rtant  effect. 

In  our  phase  screen  model  moving  the  screen  perpendicular  to  the  line  of 
sight  moves  the  intensity  pattern  with  it.  Therefore,  wind  convecting  the  tur¬ 
bulent  eddys  through  the  atmosphere  shift  the  pattern  on  the  ground.  This  idea 
is  complicated  by  the  fact  that  the  turbulent  eddys  are  themselves  rearranging 
in  time  and  that  the  atmosphere  is  not  a  phase  screen.  Both  of  these  effects  will 
reduce  the  ctarelation  between  two  intensity  detectors  on  the  ground.  When  the 
pattern  has  a  large  smal  ratio,  only  the  wind  velocity  component  perpendicular 
to  the  bands  will  be  evident  in  the  observed  motion.  If  the  wind  is  parallel  with 
the  bands,  they  may  be  very  difficult  to  see  even  if  the  intensity  fluctuations  are 
relatively  strong.  Under  such  circumstances,  the  dominant  effect  would  be  due 
to  turbulent  rearrangement.  This  would  lead  to  a  "shimmering”  effect  that 
might  not  catch  the  ejre  as  well  as  moving  bands.  Combiniiig  this  effect  with  the 
changing  orientation  of  the  crescent  could  explain  why  shadow  bands  are  some¬ 
times  seen  on  one  side  of  totality  but  not  the  other.  Annular  eclipses  and 
eclipses  with  small  values  of  e  have  much  less  of  an  anisotn^y  in  |£'(sg)|* .  In 
these  eases  the  wind  has  a  much  higher  probability  of  blowii^  in  a  direction 
parallel  to  the  crescent.  This  is  one  of  the  reasons  why  small  e  eclipses  are  not 
likely  to  generate  profound  shadow  band  events. 


The  cberecierifltic  length  scale  parallel  to  the  bands  is  roughly 
characteristic  length  scale  perpendicular  to  the  bonds 
is  not  as  straightforward  to  calculate.  1110  intensity  spectrum  perpendicular  to 
the  crescent  for  a  screen  at  a  given  altitude  has  three  characteristic  scales 
(figure  9).  The  largest  scale  is 

«2ns/?5T>/6a|(e.T}  (5. 1) 

wtuch  corresponds  to  the  low  frequency  peak  in  the  spectrum.  The  parameter 
is  given  in  eqn.  (4.4).  The  next  scale  is  where  the  source  spectrum's 
asymptotic  form  has  its  first  null. 

laMVM  =  (5-2) 

The  Uurd  scale  is  the  length  corresponding  to  the  first  nontrivial  zero  of  the 
FVesnel  filter. 

(5.3) 

Which  of  the  above  lengths  will  be  the  dominant  shadow  band  scale  will  depend 
on  the  distance  of  Uie  screen  from  the  observer  and  the  time  to  totality,  r.  The 
eye's  response  to  moving  intensity  patterns  is  a  strong  function  of  the  pattern's 
observed  angular  velocity  [lO].  Visual  sensitivity  also  drops  off  precipitously  at 
both  low  and  high  spatial  frequencies.  For  common  shadow  band  contrast  levels 
of  2  to  3%  the  lowest  visible  frequency  differs  from  the  highest  by  less  than  a  fac¬ 
tor  of  50.  bandwidth  decreasing  with  increasing  pattern  speed.  This  implies  that 
if  an  observer  can  see  scales  on  the  order  of  a  few  centimeters,  scales  larger 
than  1  to  2  meters  would  be  invisible.  Therefore,  even  thou^  there  is  a  great 
deal  more  variance  at  Lp^  than  in  the  vicinity  of  or  these  scales 

will  be  much  more  visible.  The  only  exception  is  when  the  screen  is  very  near 
the  observer,  say  100  meters  or  less.  In  such  a  case  the  Fresnel  filter  and  the 
source  spectrum  comtwe  to  form  a  high-pass  spaUal  filter  in  one  direction. 
This  will  cause  an  effect  similar  to  a  crude  Schlieren  camera  about  the  mean 
intensity.  The  effect  would  be  to  create  intensity  fiuctuations  on  the  ground  that 
are  like  refractive  "shadows”  of  the  turbulent  irregularities.  These  patterns 
would  an>oar  more  like  "smoky  wisps”  than  shadow  bands.  The  broad  source  in 


124 


the  direction  pereilei  to  the  crescent  will  limit  the  refrective  resolution  of  the 
irregulerities  to  -  Thus  the  "smoky  wisps”  would  heve  e  eherecteristic 
width  snd  a  longitudinal  scale  of  roughly  At  longitudinal  distances 

greater  than  the  pattern  would  appear  to  lose  linearity  and  wander  ran¬ 

domly.  When  the  screen  is  1000  meters  or  more  away  from  the  observer,  the 
fluctuations  with  scale  will  be  too  large  to  see  until  very  near  totality.  Hie 
scales  of  interest  are  now  end  Two  minutes  away  from  totality,  a 

screen  at  1000  meters  will  have  a  value  of  X>Mym  of  ebout  30  cm.  As  totality 
nears,  ^anve*  '*111  decrease  linearly  with  r.  A  screen  at  1000  meters  will  have  a 
Fresnel  length  scale  of  about  2.5  cm.  As  Lmwm  drops  toward  the  power  in 

the  intensity  spectrum  will  increase  in  the  vicinity  of  When  this  occurs, 

the  pattern  will  develop  structures  with  a  characteristic  scale  When 

i>  smaller  than  ^-jTrmrt  the  larger  scale  will  remain  dominant.  Whenever 
structures  at  the  Fresnel  scale  are  visible,  the  spacing  will  be  wavelength  depen¬ 
dent  since  Also,  the  intensity  variance  in  the  structure  will  be 

greater  with  smaller  wavelength.  From  the  time  when  ^  -  ^r-nr-  to  totality, 

the  pattern  will  remain  wavelength  dependent  with  scale  ^-mnir  and  the  con¬ 
trast  will  continue  to  increase.  In  the  last  few  mnnents  before  totality,  will 
drop  into  the  visible  range  and  provide  smne  larger  scale  contrast  patterns 
superimposed  on  the  smaller,  wavelength  dependent  shadow  bands.  A  screen  at 
an  altitude  of  10  km  would  have  the  same  behavior  as  given  above  except  for  the 
fact  that  source  averaging  will  keep  the  pattern’s  contrast  low  until  relatively 
near  totality.  For  a  high  altitude  screen,  it  is  only  within  a  few  seconds  of  total¬ 
ity  that  sufficient  scintillation  power  is  present  to  make  the  patterns  visible. 

How  well  the  shadow  bands  display  frequency  dependent  structures 
depends  on  e.  Ifhile  the  asymptotic  form  of  depends  only  weakly  on  e, 

the  total  intensity,  depends  strongly  on  e.  Snee  Ium  is  much  larger  for 
small  t  eclipses,  the  relative  confrosl  in  the  asymptotic  region  will  be  smaller. 
This  means  that  ecUpses  with  larger  values  of  e  will  be  better  at  mairing  the  nar¬ 
row,  high  axial-ratio  patterns  usually  associated  with  shadow  bands.  Neverthe¬ 
less,  small  e  eclipses  would  still  display  the  features  associated  with  ^ce 
fyriiv  affects  are  only  visible  when  the  screen  is  close  to  the  observer,  we  expect 
that  small  t  eclipses  or  annular  eclipses  will  predominantly  display  the  random. 
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smoky  potUrns  caused  by  low  level  Schlieren  efleeia. 

The  lest  point  to  consider  is  how  en  intensity  detector's  time  series  is 
related  to  whet  we  have  calculated.  UnfoKunately.  shadow  band  time  series  are 
very  complicated  and  require  a  great  deal  of  information  in  order  to  analyze. 
The  simplest  approach  is  to  assume  that  the  turbulence  is  "frozen”  and  does  not 
rearrange  as  the  wind  moves  it  past  (Taylor's  hjrpothesis).  As  such,  the  intensity 
pattern  moves  at  a  velocity  equal  to  the  wind  component  perpendicular  to  the 
line  of  sight.  For  the  moment,  we  will  ignore  the  fact  that  the  source  is  changing 
rapidly  in  time.  The  wind  velocity  vector  perpendicular  to  the  line  of  sight  is  f. 
Using  Taylor's  hypothesis  we  see  that 

/(t<o)  =  /(t*ff.fo  +  0 

which  leads  us  to  the  temporal  intensity  correlation  as 

-  Q(Tf)  (5.4) 

lyansforming  (5.4)  with  respect  to  f  to  find  the  temporal  intensity  power  spec¬ 
trum  we  see  that 

(S.5) 

where  cg  and  are  the  components  of  d  parallel  and  perpendicular  to  the  wind 
vector,  T,  and  /  is  the  temporal  frequency  in  cycles  per  second.  We  may  include 
the  effect  of  the  finite  detector  aperture  by  multipl3nng  the  right  band  side  of 
eqn.  (2.27)  by  (2n)*|7(2)f  idiere  X{i)  is  the  Fourier  transform  of  the  detector's 
aperture  distributim.  The  difficulty  in  the  analysis  comes  b-om  the  fact  that  the 
wind  vector  typically  lies  alrnig  some  direction  other  than  the  perpendicular  to 
the  crescent.  This  leads  to  a  rather  complicated  integral.  The  "strip  integra¬ 
tion"  in  (5.5)  washes  out  most  of  the  oscillatory  structure  in  the  source  spec¬ 
trum  and  the  Fresnel  filter  allowing  us  to  find  approziniate  scafing  laws  for  eipi. 
(5.5)  when  the  wind  is  perpendicular  to  the  shadow  bands.  For  temporal  frequen¬ 
cies  greater  than  *•  ^  ^nnriMn  where  is  the  size  of  the  detector  ^>erture, 

there  is  essentially  no  power.  If  "  '  snd  the  apparent 

wind  is  at  angle  less  than  ^‘9s  ^  axis,  the  temporal  intensity 


tpectrum  will  ical«  approiiiMUly  m  In  tbt  frtqiMney  range 

^  icalini  ^  roughly  For 

y/ <liv/ Lmmvt  or  *tion  Iawmi  >  ^wivw  the  aptictrum  will  depend 
strongly  on  the  details  of  the  source  spectrum  and  the  wind  direction.  For  very 
low  temporal  frequencies,  f<cu/Iyj^.  the  spectral  power  will  scale  approxi¬ 
mately  as  If  the  apparent  wind  is  at  an  angle  greater  than  ^'fs  the 

lower  temporal  frequencies  are  greatly  enhanced.  For  this  case,  the  analysis  of 
the  temporal  spectrum’s  shape  would  require  detailed  numerical  calculations 
specific  to  the  observation  geometry  and  the  wind  direction. 

The  inherent  difficulties  in  shadow  band  time  series  analysis  are  com¬ 
pounded  by  the  dynamic  character  of  the  source  distribution.  Estimates  of  the 
tempwal  intensity  spectrum  must  be  done  on  time  series  whose  length  is 
sufficiently  short  that  the  source  distribution  may  be  considered  unchanging. 
IVom  this  perspective  it  may  be  preferred  to  examine  time  series  taken  rela¬ 
tively  far  from  totality.  Unfortunately,  this  does  not  permit  making  use  of  the 
stronger  fluctuaticms  near  totality.  The  other  qption  would  be  to  use  shorter 
time  series.  However,  this  approach  would  lead  to  prohibitive  spectral  estima¬ 
tion  error  unless  the  wind  speed  was  quite  high.  These  problems  are  particu¬ 
larly  difficult  if  we  try  to  observe  the  low-frequency  spectral  behavior  for 

8.  Use  BBecta  of  Extended  Atmoapherlc  Turbulence. 

Atmospheric  turbulence  is  not  confined  to  a  thin  lajrer.  It  is  (Nresent  every¬ 
where  from  the  ground  to  tens  oS  kiiometers  in  altitude.  The  models  described 
in  the  previous  sections  need  to  be  modified  to  take  this  into  account  As  long 
as  the  scattering  is  weak,  the  intensity  fluctuatim  spectrum  is  the  linear  super¬ 
position  of  the  thin  screen  results.  It  is  given  by 

♦/•(<)  -  (2n)^irk*7*,»(*;s)|J^(sg)|*sin*(^^(is 
0  ** 

vdiere  is  the  turbulence  spectrum  at  a  distance  s  from  the  observer.  If 

the  sun  is  not  at  the  xenith,  the  altitude  is  given  by  s  cos(0  where  is  the  sun's 
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senith  uifle.  The  reiulting  formula  for  the  intamity  tpectrum  bceomat 

♦/(*)  *  ♦  {2«mi38)«  *7 qf(«coe(C))ein*(^d«  (6.1) 

where  ^(A)  is  the  turbulence  structure  cmiatant  at  altitude  h.  Combining  eqs. 
(6.1),  (4.5),  and  (2.17)  gives  the  total  scintillation  index  as 

m*(X,e,T,0  =  J^CSiMcosiOyfl^ii*.Ke.T)d*  (6.2) 

Hie  integrand  in  (6.2)  shows  the  relative  importance  of  each  distance  to  the 
total  intensity  fluctuation  spectrum.  Given  the  turbulence  strei^th  proflle, 
CJf(h),  we  can  see  which  altitudes  are  contributing  meet  to  the  shadow  bands. 
Cli(h)  is  known  to  vary  with  time  of  day.  season,  geographical  position,  humidity, 
wind  speed,  and  surface  texture  [11].  For  our  purposes,  we  will  use  a  comlana* 
tion  of  models  to  generate  a  generic  turbulence  profile.  Our  proflle  will  be  a 
form  of  the  Hufhagel  model  [12],  modified  to  account  for  the  strong  turbulence 
at  low  altitudes,  given  by 

q*(/k)  s  q?(l)/i-^»e-*^'«>  +  (^7xl0-*•)e (6.3) 

where  h  is  the  altitude  in  meters,  C«(l)  is  the  turbulence  structure  constant  at 
1  meter,  and  <i^>  is  the  mean-square  wind  speed  at  high  altitudes.  For  the  cal¬ 
culations  that  follow,  we  asrune  ''typical”  values  for  the  constants  in  (6.3): 
Cm(1)  *  6xi0“'*m~*'*  and  <i/*>  «  850bi*/s*.  The  resulting  turbulence  proflle  is 
shown  in  figure  10.  The  peak  at  lOOOO  m  is  the  enhanced  turbulence  at  the  tro- 
popause.  Ihe  resulting  contribution  to  the  total  scintillation  index  for  various 
altitudes  is  shown  in  figure  11.  Figure  11a  shows  the  temporal  evolution  of  the 
m*  proflle  for  a  t]fpical  eclipse.  At  about  two  minutes  before  totality,  only  tur¬ 
bulence  very  near  the  ground  is  making  substantial  contributions  to  the  inten¬ 
sity  fluctuations.  Even  by  20  seconds  before  totality,  the  turbulence  within  the 
bottom  two  kilometers  is  responsible  for  essentially  all  of  the  scintillations.  Only 
within  5  to  10  seconds  of  totality  does  the  high  altitude  turbulence  make  a 
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noUcMbIc  tllect.  Hoirevtr,  •v«n  at  this  late  time,  the  turbulence  within  two 
kilometers  of  the  ground  is  much  more  important  than  the  turbulence  at  the 
tropopause.  Figure  11b  shows  the  dependence  of  the  scintillation  profile  on 
wavelength.  At  twenty  seconds  before  totality,  the  wavelength  dependence  is 
confined  to  within  a  few  hundred  meters  of  the  ground.  By  ten  seconds  away 
from  totality,  the  wavelength  dependence  has  spread  to  include  all  of  the  lower 
altitude  turbulence  while  the  scintillations  from  the  tropopause  are  still 
wavelength  independent.  Only  within  a  few  seconds  of  totality  does  the  tropo¬ 
pause  produce  wavelength  dependent  scintillations.  Figure  11c  shows  that  the 
scintillation  contribution  is  degraded  at  all  altitudes  equally  for  smaller  e 
eclipses.  Figure  lid  shows  the  effect  of  zenith  angle  on  the  scintillation  profile. 
Note  that  not  only  are  the  scintillations  stronger,  they  are  the  result  of  more 
distant  turbulence.  The  total  scintillation  index  is  shown  in  figure  3.  An 
observer's  estimate  of  the  intensity  variance  would  depend  on  the  scales 
observed.  The  calculations  shown  in  figure  12  include  oil  scales.  The  wavelength 
dependence  shown  in  figure  12a  would  be  much  more  pronounced  if  the  large- 
scale.  wavelength  independent  scintillations  were  ignored.  The  same  point  is 
true  of  figure  12b  which  would  show  a  much  greater  decrease  in  for  smaller 
values  of  e  if  the  large-scale  irregularities  were  1^  out.  The  dependence  vl  m* 
on  zenith  angle  is  shown  In  figure  13.  As  long  as  the  turbulence  in  the  bottom 
few  kilometers  dominates  the  scintillations,  we  may  treat  CS  as  being  propor¬ 
tional  to  A~^^.  This  leads  to  an  af^roximate  scaling  formula  for  m*  with  zenith 
angle  as 

m*(0  =  (.ec(Or*m«(C  =  0)  (6.4) 

The  m*  profile  allows  us  to  define  a  ’'tjrpical”  scattering  distance  which  is 
the  centroid  of  the  scintillation  contributions.  We  define  the  n-tb  moment  of  the 
m’  profile  as 


/s"  qf(s  coo(f  ))q.^s  )ds 

5 _ 

^qf(scoa(0)n»s(«)^ 


(6.5) 


129 


« 


W 


Th«  txptcal  icatUrinc  distance  ii  f  and  the  typical  width  of  the  scatterinf 
region  la 

C  = 

A  more  meaningful  estimate  of  the  typical  scattering  distances  is  achieved  by 
treatii^  the  low  altitude  and  the  tropopause  turbulence  contributions 
separately.  In  our  turbulence  proflle  model  eq.  (6.3),  the  tropopause  tur¬ 
bulence  is  described  by  the  last  term.  Using  the  moments  of  these  terms 
separately,  two  typical  scattering  distances  are  found,  figure  14.  Note  that  the 
typical  scattering  distance  of  the  low  altitude  turbulence  is  rather  insensitive  to 
zenith  angle  while  the  tr(^>opause  distance  increases  with  C. 

The  structure  of  shadow  bands  in  the  extended  turbulence  model  is  a  com- 
l^ation  of  thin  screen  model  structures  from  all  of  the  contributing  distances. 
A  good  picture  of  the  evolving  scales  may  be  found,  without  making  a  detailed 
spectrum  calculation,  by  using  the  thin  screen  model  appropriate  for  the  typical 
scattering  distance.  The  structure  is  found  by  using  eqs.  (5.1),  (5.2).  and  (5.3) 
with  s  replaced  by  i.  The  effect  of  the  actual,  distributed  turbulence  will  be  to 
make  the  thin  screen  scales  less  ihstinct.  rather  than  to  introduce  any  new 
scales.  The  thin  screen  scales  are  shown  in  figure  15.  Features  much  larger 
than  1  meter  are  not  visible  and  are  not  shown  although  they  contribute  to  m^. 
When  the  scintillations  first  become  visible,  they  present  wavelength  indepen¬ 
dent  structures  with  scale  frmn  turbulence  very  near  the  ground.  As  total¬ 
ity  nears,  the  tjrpical  scattering  distance  goes  up  and  the  prominent  scale 
becomes  Lmmtu  lo^  altitude  turbulence.  The  scintillations  from  the 

tropopause  are  still  suffering  from  source  averaging  and  are  too  faint  to  be  seen. 
As  the  crescent  narrows,  the  axial  ratio  of  the  pattern  increases  and  the  pat¬ 
terns  become  true  "shadow  bands.”  The  very  low  pattern  contrasts  will  lead  to 
the  observation  that  the  bands  are  narrow  dark  strips  separated  by  larger 
bright  regions  due  to  the  logarithmic  sensitivity  of  the  eyes  [13].  As  the  con¬ 
trast  increases,  the  pattern  will  appear  to  have  roughly  equal  strips  of  light  and 
dark.  Finally,  the  structures  at  scale  merge  with  stractxires  at  scale 

/■Pi-  and  become  wavelength  dependent  Once  the  dominant  scale, 

the  pattern  spacing  will  slowly  increase  as  the  t]rpical  scattering  distance 


continual  to  incrioso.  Moinwhile.  icintiUoUon  itructuros  with  icilo  will 
bov«  returned  to  a  vtiible  aize,  diaptajring  email  ajdal*ratio.  wavelenf  th  indepen* 
dent  patterns  superimposed  on  the  shadow  bands.  In  the  last  few  seconds  from 
totality,  the  tropopause  finally  makes  a  visible  effect,  with  patterns  at  scale 

rapidly  merging  with  the  tropopause's  f  ~ -  If  the  eclipse  has  a  small 

value  of  e,  the  axial  ratio  is  never  as  great  as  in  longer  eclipses.  In  additioa  the 
mean  intensity  is  greater  causing  a  diminished  contrast  at  all  scales  except 
Lp^.  Therefore,  small  t  eclipses  and  annular  eclipses  will  predominantly  gen¬ 
erate  the  randmn.  smoky  patterns  associated  with  Figwe  15b  shows  the 
effect  of  zenith  angle  on  the  characteristic  scales.  While  the  low  altitude  struc¬ 
tures  are  relatively  unaffected,  the  turbulence  at  the  tropopause  remains 
heavily  source  averaged  until  even  closer  to  totality. 

Since  the  wind  velocity  is  often  a  function  of  altitude,  the  shadow  band 
motion  may  change  as  the  typical  scattering  distance  moves  to  greater  alti¬ 
tudes.  If  a  wind-shear  layer  is  present,  distinct  shadow  band  patterns  may  be 
seen  superimposed  on  each  other,  moving  with  different  speeds.  If  the  observer 
is  not  on  the  center-line  of  the  path  of  totality,  the  band  motion  wiU  be  even 
more  complex:  the  wind  velocity  component  perpendicular  to  the  bands  chang¬ 
ing  as  the  crescent  rotates.  The  complications  caused  by  wind  speed  variations 
are  jret  another  insurmountable  problem  in  shadow-band  time  series  analjfsis. 
The  most  reasonable  ^proach  would  be  to  treat  the  scinUUatioos  as  if  they 
arose  from  a  thin  screen  located  at  the  typical  scattering  distance  and  ignore 
changes  in  the  shadow  band  speed  during  the  observation.  Nevertheless,  the 
error  analysis  should  seriously  take  such  factors  into  account 

The  generic  turbulence  profile,  eq.  (6.3).  has  given  us  a  perspective  on  the 
evolution  of  a  typical  shadow  band  event  However,  the  actual  turbulence  profile 
present  at  an  observation  site  may  (hfler  greatly  from  the  modeL  In  adthtion  to 
large-scale  turbulence  distributions,  an  actual  C!  profile  is  nof  smooth;  it  is  lay^ 
ered  and  patchy[14].  One  general  principle  may  be  inferred  from  the  generic 
prt^e:  shadow  band  visibility  is  largely  controlled  by  the  level  of  turbulence 
within  the  bottom  2-3  kilometers  of  the  atmosphere.  Therefore,  any  situation  in 
which  the  lower  altitude  turbulence  is  diminished  will  not  present  a  strong  sha¬ 
dow  band  display.  Such  cases  would  include  sites  with  very  good  "seeing"  or 
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plAccs  with  high  humidity,  such  sm  near  lakes  or  on  the  ocean.  Observations 
made  under  very  windy  conditions  shotild  not  make  good  visual  shadow  band 
events  due  to  the  response  time  of  the  eye.  However,  properly  sampled  pho* 
toelectric  observations  obtained  under  similar  circumstances  and  would  lead  to 
better  temporal  spectrum  estimates.  Sites  where  the  eclipse  is  seen  further 
from  the  zenith  will  show  stronger  shadow  bands. 


In  the  preceding  sections,  we  have  seen  how  weak  scattering  scintillation 
theory  can  be  used  to  explain  all  of  the  salient  features  of  eclipse  shadow  bands. 
The  unknown  turbulence  and  wind  speed  profiles  are  the  only  reasons  irtiy  expli¬ 
cit  shadow  band  predictions  are  not  possible.  However,  with  reasonable  assump¬ 
tions  about  the  meteorology,  many  features  of  the  shadow  bands  may  be  under¬ 
stood.  The  band  motion  is  driven  by  the  wind  velocity  component  at  various  con¬ 
tributing  altitudes  perpendicular  to  the  bands.  The  band  orientation  is  parallel 
with  the  crescent  which  is  a  function  of  time  if  the  observer  is  not  on  the  center 
of  the  path  of  totality.  The  bands  will  be  perpendicular  to  the  lunar  shadow's 
path  far  from  totality  and  parallel  with  edge  of  the  shadow  near  totality.  The 
band  spacing  is  a  complicated  interplay  erf  three  intensity  spectrum  scales  and  a 
changing  scattering  altitude.  The  patterns  staK  out  random  and  disorganized, 
becoming  more  linear  and  organized  as  tbe  crescent  narrows  and  the  scattering 
height  increases.  The  band  spacing  decreases  as  totality  nears  and  finally 
becomes  wavelength  dependent  with  a  spacing  proportional  to  VX.  The  pattern 
contrast  is  wavelength  dependent  as  diffraction  effects  become  more  in^xxlant 
near  totality,  shorter  wavelengths  exhibiting  higher  cmitrast  patterns.  Both  pat¬ 
tern  contrast  and  spacing  increase  with  increasing  zenith  angle.  The  scintilla¬ 
tion  theory  provides  no  means  for  polarization  effects  and  none  are  expected. 
Shorter  eclipses  and  annular  eclipses  are  not  effiment  at  generating  high  axial- 
ratio  shadow  bands:  only  random,  smoky  patterns.  Shadow  bands  are  related  to 
the  same  turbulence  responsible  for  "seeing":  good  seeing  in4>lies  poorer  sha¬ 
dow  band  contrast.  All  shadow  band  phenomena  are  statistically  symmetric 
before  and  after  totality,  except  for  motion.  The  temporal  structure  of  the 
intensity  fluctuaUons  is  strongly  dependent  on  the  direction  and  altitude 
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distribution  of  tho  wind. 

Sevorol  txperimontal  ttudiea  at  tho  thodow  bond  pbenomono  hovo  boon 
porformod.  Unfortunatoly,  publiihod  oxporimonU  with  sufficiont  information 
and  quality  to  mako  usoful  theoretical  comparisons  are  rare.  Three  of  the  best 
are  Hults  tt  al  [15].  Quann  and  Daly  [4],  and  Marschall  tf  al  [16].  Hults  sf  ai 
made  a  study  of  the  band  orientation  within  the  path  of  totality.  They  note  a 
"curious  rotation"  of  the  bands  when  the  observer  was  not  on  the  center^line  of 
the  path  of  totality.  Quann  and  Daly  noted  the  same  rotation  but  correctly  attri¬ 
buted  it  to  the  apparent  rotation  of  the  solar  crescent  Their  position  was  just 
outside  of  the  path  erf  totality  and  hence  would  have  observed  the  solar  crescent 
rotate  by  nearly  180*.  They  reported  data  from  a  rawinsonde.  launched  prior  to 
the  eclipse,  which  showed  the  wind  to  be  roughly  parallel  with  the  path  of  total¬ 
ity.  They  also  made  a  detailed  photoelectric  survey  to  study  polarization  and 
wavelength  effects.  They  found  no  p<rfarizatioo  effects  and  their  aggregate 
power  spectra  have  widths  which  scale  roughly  as  VX.  Unfortunately,  they 
apparently  averaged  all  times  into  their  spectra  and  no  real  conclusions  may  be 
drawn.  In  the  blue  and  ultraviolet,  they  observed  the  bands  to  slow  down  and 
become  invisible  near  maximum  partial  eclipse.  This  was  not  seen  in  the  green- 
yellow  observations.  During  the  Ume  of  this  "cutcrfl,"  the  wind  was  roughly  paral¬ 
lel  to  the  bands  and  the  wavelengths  affected  would  have  had  the  greatest 
ratia  Additionally,  the  actual  effect  of  solar  limb  darkening  would  be  to  reduce 
the  mean  intensity  of  the  solar  crescent  in  blue  and  ultraviolet  near  totality. 
Therefore,  the  observed  cutoff  may  have  been  due  to  insufficient  mean  intensity 
at  those  wavelengths.  The  best  intensity  fluctuation  spectrum  measurements  to 
date  were  made  by  Marschall,  Mahon,  and  Henry  [16].  Their  two  published  spec¬ 
tra  were  taken  sufficiently  far  from  tcrfality  that  the  source  evolution  could  be 
ignored.  The  noise  in  the  pre-totality  spectrum  appears  more  pronounced  than 
in  the  post-totality  spectrum,  suggesting  that  the  wind  was  more  nearly  perpen¬ 
dicular  to  the  bands  after  totality.  Assuming  v/Ln^g^u  to  be  about  30  Hz,  the 
approximate  scaling  form  of  the  power  spectrum  may  be  flt  to  the  post-totality 
spectrum  (flgure  16).  The  thin  screen  spectrum  is  too  steep  above  v/^Awmi* 
This  is  doubtless  due  to  the  distributed  nature  of  the  turbulence,  lower  altitudes 
providing  higher  frequency  intensity  fluctuations. 


Despite  the  high  quality  of  the  power  speetnim  obtained  by  MarschaL  at  at, 
the  uncertainties  involved  in  the  time  dependence  and  the  wind  direction  make 
an  exact  comparison  with  the  theory  dUTicult  Two  dimensional  observations 
that  did  not  depend  on  pattern  drift  would  be  far  better.  Perhaps  with  the 
advent  of  higher  quality  video  cameras  and  photographic  film,  better  comparis¬ 
ons  will  be  possible  in  the  future. 
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Flfun  CapUooa 

1.  Geometry  for  the  thin  screen  model. 

2.  Evolution  of  solar  crescents  for  an  observer  on  the  center  of  the  path  of 
totality.  Crescents  are  shown  at  r  intervals  of  0.01  and  e  is  (a)  0.01,  (b) 
0.06,  and  (c)  0.10. 

3.  Definitions  of  crescent  angles  and  lengths. 

4.  Evolution  of  crescent  horn  angles  for  t  -  0.001,  0.01,  0.06,  and  0. 1. 

5.  Mean  intensity  of  solar  crescent  for  e  =  0.001,  0.01,  0.05,  and  0.1.  Intensity 
has  been  normalized  to  the  full  disk  intensity. 

6.  Evolution  of  solar  crescents  for  an  observer  not  on  the  center  of  the  path  of 
totality.  Crescents  are  shown  at  T  intervals  of  0.01  for  e  =  0.06  and  (a) 
p  =  0.9,  and  (b)  p  =  0.6. 

7.  An  example  demonstrating  the  relationship  between  the  actual,  normalized 
time  to  totality,  Tq.  and  the  geometrical  t.  The  case  shown  is  for  it  =  0.05 
and  p  =  0  to  1.2  at  intervals  of  0. 1. 

8.  A  slice  through  the  Fourier  transform  of  the  crescent  along  the  axis  for 
(a)  e  +  T  =  0.1.  (b)  c  +  t  =  0.01.  and(c)  e  +t= 0.001. 

9.  Scintillation  efficiency  profiles  vs.  distance  from  the  observer.  Figure  9a 
shows  the  evolution  of  the  efficiency  profile  vs.  r.  Figure  9b  shows  the  effect 
of  changing  the  wavelength  for  three  different  times.  Figure  9c  shows  the 
efficiency  profile  for  three  different  values  of  e  at  two  different  times,  t. 

10.  Plot  of  the  shadow  band  intensity  fluctuation  spectrum  perpendicular  to  the 
solar  crescent  showing  the  three  major  scales. 

11.  Plot  of  the  modified  Huhiagel  turbulence  profile. 

12.  Scintillation  index  profile  evolution. 

13.  Scintillation  index  evolution. 

14.  Dependence  of  the  scintillaUon  mdex  on  zenith  angle. 

15.  Typical  scattering  distances  vs.  time  and  zenith  angle. 

16.  Evolution  of  the  characteristic  scales  vs.  time  and  zenith  angle. 

17.  Observed  shadow  band  time  series  power  spectrum  compared  to  the 
approximate  thin  screen  model 
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